Double-component convection due to different 
boundary conditions in an infinite slot 
diversely oriented to the gravity 1 

N. Tsitverblit 2 

Department of Fluid Mechanics and Heat Transfer, Tel-Aviv University, 

Ramat-Aviv 69978, Israel 

Abstract 

Onset of small-amplitude oscillatory and both small- and finite-amplitude steady 
double-component convection arising due to component different boundary con- 
ditions in an infinite slot is studied for various slot orientations to the gravity. 
The main focus is on two compensating background gradients of the components. 
The physical mechanisms underlying steady and oscillatory convection are analyzed 
from the perspective of a universally consistent understanding of the effects of dif- 
ferent boundary conditions. In a horizontal slot with inviscid fluid addressed by 
Welander [Tellus, Ser. A 41 (1989) 66], oscillatory convection sets in with the most 
unstable wave number and oscillation frequency being zero. Exact expressions for 
the critical fixed-value background gradient and the respective group velocity at 
zero wave number are derived from the long-wavelength expansion both for the 
horizontal slot with independently varying background gradients and for the in- 
clined slot with the compensating gradients. In the horizontal slot with viscous 
fluid, the dissipation of along-slot perturbation-cell motion reduces efficiency of the 
oscillatory instability feedback and thus prevents the most unstable wavelength 
from being infinite. Based on this interpretation, the oscillatory instability of a 
three-dimensional (3D) nature is predicted for an interval of long two-dimensional 
(2D) wavelengths in an inclined slot, and such 3D instability is indeed shown to 
arise. Related general conditions for three-dimensionality of most unstable distur- 
bances are also formulated. As the slot orientation changes from the horizontal 
by angle 6 ( > vr/2), the oscillatory 2D marginal-stability boundaries in invis- 
cid and viscous fluid are expected to eventually transform into respective steady 
ones. Oscillatory instability in the vertical slot with viscous fluid, first reported 
by Tsitverblit [Phys. Rev. E 62 (2000) R7591], is of a quasisteady nature. Its (new) 
mechanism is identified. It is underlain by differential gradient diffusion. As the hori- 
zontal slot at 9 = 7r, addressed by Tsitverblit [Phys. Fluids 9 (1997) 2458], changes 
its orientation to vertical, the wave number interval of linear steady instability 
shrinks to the vicinity of the most unstable zero wave number and vanishes. Con- 
sistently with the basic nature of finite-amplitude steady convection being the same 
in the horizontal and vertical slots, the respective convective flows are continu- 
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ously transformed into each other. The dissimilarity between the nature of finite- 
amplitude steady convective flows in the horizontal slot with 6 = 0, revealed by 
Tsitverblit [Phys. Lett. A 329 (2004) 445], and that in the vertical slot is shown to 
eventually give rise to a region of hysteresis in 6 £ (0, n/2). 
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1 Introduction 



This work, some of whose aspects were promulgated in pQ , addresses onset of 
double-component, buoyancy-driven convection resulting from the boundary 
conditions for one component being different from those for the other. Such 
convection has recently been identified as a fundamentally new class of hydro- 
dynamic instabilities underlying the formation of spatial and temporal flow 
patterns from a steady equilibrium state of spatially homogeneous fluid. The 
main emphasis of this study is on developing a universally consistent perspec- 
tive for understanding the effects of different orientation of the fluid domain 
to the gravity. 

Double-component convection is relevant to the phenomena in small-scale 
oceanography [2], geology [3], geodynamo [I], and such areas of astrophysics 
as ordinary evolution of stars [5] and core-collapse supernova explosions [6]. 
It has technological applications as well, in particular to crystal growth [TJ. 
Double-component convection has also been suggested as a possible cause of 
layering phenomena in colloidal suspensions [S] and as the origin of pattern 
formation in soap films [9] . Early reviews of some of the above applications of 
the effects of different component diffusivities in double-component convection 
can be found in [TO] . 

In addition, convective flows are commonly used in fundamental studies of 

1 NOTICE: this is the author's version of a work that was accepted for publication 
in Annals of Physics. Changes resulting from the publishing process, such as peer 
review, editing, corrections, structural formatting, and other quality control mecha- 
nisms may not be reflected in this document. Changes may have been made to this 
work since it was submitted for publication. A definitive version was subsequently 
published in Annals of Physics [Ann. Phys. (2007) 322(8) 1727-1770] 

DOI 10.1016/j.aop.2006.10.001. 

2 Address for correspondence: 1 Yanosh Korchak Street, apt. 6, Netanya 42495, 
Israel; e-mail: naftali@eng.tau.ac.il 



2 



transition to turbulence and nonlinear pattern formation [TTlfT2] . Double- 
component flows where a distinction between the components comes from 
component different boundary conditions are also of basic significance in the 
context of their applications to large-scale environmental and turbulent pro- 
cesses. In particular, such processes are Langmuir circulations [13J and the 
global ocean thermohaline circulation [T^ffT5] . Except for limited aspects of 
the latter context, however, two-component convection in pure fluid has been 
previously addressed mostly in the framework of only the effects of different 
diffusion coefficients [TBI 1711 181. 



Welander [19] was the first to highlight the onset of oscillatory convection re- 
sulting from an unequal effect of different boundary conditions on the compo- 
nent diffusion gradients in perturbed state. Such convection was demonstrated 
in [Tj5] to arise in a horizontal layer of inviscid fluid with a double-component, 
statically stable net stratification. The instability mechanism described in [TjJ] 
has later been recognized by the present author [201121] as being conceptually 
analogous to that driving oscillatory convection in the diffusive regime of the 
classical double-diffusion [T6]fT7] . It has also been suggested in [20ll21f22j that 
such analogy between the effect of different boundary conditions and that of 
unequal diffusivities is of a generic nature. Near a boundary towards which 
across-slot perturbation motion is directed, the component with its bound- 
ary value fixed would have a higher perturbation gradient than the compo- 
nent whose boundary condition is specified in terms of the flux. The different 
rates of gradient diffusion resulting from such a disparity could thus trigger 
double-component convection analogously to the effect of unequal component 
diffusivities in conventional double-diffusive convection. 

In terms of the analogy just described, the mechanism of steady convection un- 
der the stratification inverse to that in [Tj5] has been shown in [2"U1l2~Tj to be con- 
ceptually reminiscent of the finger instability in the classical double-diffusion 
[TB] . (The mathematical problem discussed in [20f2TJ also happens to describe 
double-component Langmuir circulations [13].) In application to two horizon- 
tal gradients, the effect of different boundary conditions was also demonstrated 
in [22] to apply to a laterally heated stably stratified slot (LHSSS). The lat- 
ter configuration had been studied before only in the context of the classical 
double-diffusion [2"3] . 

The above effect of boundary conditions was also demonstrated in |24j to result 
in steady convection when two compensating horizontal gradients maintained 
by different boundary conditions are applied to a vertical layer of Boussinesq 
fluid at rest. In the latter configuration, however, finite- amplitude convective 
steady flows arise without the respective linear instability of the conduction 
state, as in the scenario first proposed in |25j. This is qualitatively different 
from the analogous problem of the classical double-diffusion [26], where the 
respective linear instability is concomitant of subcritical steady convection. 
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The finite-amplitude manifestation of convection uncovered in [23] was also 
found to take place under a stable vertical solute stratification. In addition, 
the scenario in [23] unveiled the existence of an oscillatory linear instability. It 
thus suggested the possibility of a new oscillatory manifestation of convection 
due to different boundary conditions. 

Differential diffusion due to unequal perturbation gradients of the compo- 
nents (differential gradient diffusion [19.20.21,22,24]) is not the only process 
by means of which different boundary conditions can trigger convection from 
the state of rest. As reported in [TP?] , finite-amplitude steady convection 
could arise well before the onset of the respective linear instability in the 
viscous version of the problem considered in [19J. In this case, the disparity 
between component stratifications resulting from finite-amplitude perturba- 
tion generates convective flows due to the feedback coming from nonlinear 
Rayleigh — Benard convection. Such finite-amplitude convection was found in 
[2T] to coexist with an oscillatory linear instability in viscous fluid. 

The main objective of the present study is to provide a comprehensive under- 
standing of the manifestation of infinite-slot convection resulting from different 
boundary conditions for diverse slot orientations to the gravity. The central 
issue in achieving this objective is oscillatory small-amplitude manifestation 
of such convection. Resulting from this work, the understanding of this issue 
has required uncovery and analysis of inviscid and viscous transformations of 
the oscillatory instability into respective steady one. In addition, it required 
identification of a new physical mechanism of oscillatory instability. Such insta- 
bility arises in the vertical slot with viscous fluid [23]- The key physical effects 
underlying transformation between the inviscid mechanism of horizontal-slot 
oscillatory instability in [TjJ] and the new vertical-slot mechanism in viscous 
fluid are thus also described in this work. Implications of the provided inter- 
pretations, in particular for the existence of a three-dimensional oscillatory 
instability in an inclined slot, are analyzed as well. This study also completes 
the understanding of steady convection due to different boundary conditions 
when a horizontal slot is transformed into an inclined and vertical slot. 



2 The problem formulation and solution procedures 

2.1 The problem and governing equations 

Let compensating gradients of the contributions of two components to the den- 
sity be maintained by different boundary conditions at the walls of a slot with 
pure fluid whose orientation to the gravity is characterized by angle 9. This is 
a particular case of the problem illustrated in Fig. [TJ where 9 ( > in Fig. [T]) is 
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the angle between the direction opposite to the gravity and that of the across- 
slot coordinate axis. The component gradients in Fig. [1] are represented by the 
Rayleigh numbers Ra = gaATd 3 / kv and Ra s = —gf3(dS/dx)d 4 /Kv = fiRa. 
Here, x is the (dimensional) across-slot coordinate, d is the width of the slot, 
AT is the (dimensional) difference between the values of temperature (stand- 
ing for the component with fixed-value boundary conditions) at the bound- 
aries with smaller and larger across-slot coordinates, dS/dx is the boundaries- 
prescribed (dimensional) derivative of solute concentration (standing for the 
component with flux boundary conditions), a is the coefficient of thermal ex- 
pansion, (3 is the coefficient of the density variation due to the variation of 
solute concentration, g is the gravitational acceleration, v is the kinematic vis- 
cosity, and k = K7 1 = Kg is the diffusivity of both components. The bar means 
that the respective variable is dimensional. Unless explicitly stated otherwise, 
Ra > and Ra s > as well as // = 1 are assumed. 

The component diffusivities are set equal to eliminate the effects of the clas- 
sical double-diffusion. The effects of different boundary conditions could thus 
be analyzed separately from those of unequal diffusivities. Such approach is 
analogous to that used in previous studies of conventional double-diffusive con- 
vection. In most such studies, the components with unequal diffusivities have 
not been distinguished from each other in terms of boundary conditions. As in- 
dicated in [28], equal diffusivities can in principle be experimentally modeled 
with two solutes. The Prandtl number would then be significantly different 
from the Pr = 6.7 specified in the caption of Fig. [TJ This parameter, how- 
ever, is not expected to have a qualitative effect on the results and physical 
interpretations reported herein. Equal diffusivities could also be interpreted as 
eddy transport coefficients, as in [T3][T5] . In this case, Pr = 6.7 is not outside 
the range of realistic values corresponding to all diffusion coefficients being of 
the eddy type. 

The condition that the component gradients be exactly compensating (/i = 1) 
eliminates the along-slot motion arising when the slot orientation differs from 
horizontal. (This condition has also been adopted in many previous studies 
of conventional double- diffusive convection [2~6f2l?| .) In the absence of such 
a background flow, comparison of the results for an inclined or vertical slot 
with those for the horizontal slot is substantially facilitated. The same could 
also be said about the comparison between the viscous and inviscid problems, 
for the along-slot base flows would have been dissimilar in such problems. In 
particular, these comparisons are relevant for understanding the physics of 
transformation between oscillatory convection in inviscid fluid at 9 = [T5] 
and that in viscous fluid at 6 = 71 /2 |24j . 

The equations describing the two-dimensional (2D) problem in Fig. [T]in which 
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an along-slot solute stratification is also present can be written as follows: 



d( dip d( dip dC, 1 dt ds . 
dr dx dy dy dx Pr dx dx 

+ ^ + ^ (1) 
Pr dy dy dx 2 dy 2 ' 

^ = ^± + dhp_ 

dx 2 dy 2 ' 

dt_ (fy(H_(fydt_ J_ ( &t_ dHs 

dr dx dy dy dx Pr dx 2 dy 2 ' 

ds ^ dip ^ds r A ^ ® s * (® 2 ' s 4- ^ 2,s ) (4) 

dr dx dy dy dx Pr dx 2 dy 2 ' 



where the across-slot, u, and along-slot, v, velocities are 



u 



vorticity 

solute concentration 
temperature 



dip dip 
dy 1 dx' 



dv du 
dx dy ' 



S = —Ra s y + s, 



1 2 + ' 



Pr = v/k is the Prandtl number, r is the time, x G (—1/2,1/2), y G 
(— A/2, A/2), and A = X/d is the specified along-slot period. 

Unless explicitly emphasized otherwise, the along-slot solute stratification, 
characterized by the Rayleigh number Ras = gP\dS/dy\d A / ku, is assumed 
to be zero. The formulation with Ras ^ is used only for facilitating the 
discussion related to results in [22f2i] . where the problems with Ras ^ have 
been addressed. 

The above equations were considered along with the no-slip formulation of 
wall boundary conditions 

j. d 2 ip , Ra 
ds 

— = -fjtRa = -Ra s (x = =f1/2, -A/2 < y < A/2) (5) 
dx 
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and periodic boundary conditions in the along-slot direction 



£(x,A/2)=£(z,-A/2) 




(-1/2 < x < 1/2), 



s(x,X/2) = s(x,-X/2) 



ds(x,X/2) 
dy 



ds(x,-X/2) 
dy 



(-l/2<x<0, 0<x<l/2), s(0,±A/2)=0, 



(6) 



where £ stands for (, ip, and t. As indicated in the caption of Fig. [T], Pr = 6.7 
was used throughout the present study. As at the middle points of across-slot 
boundaries in (jBJ), specification of the values of s is needed to identify the 
solute concentration scale and select the phase of a nontrivial flow. Technical 
difficulties that still exist with the employed approach will be discussed in Sec. 



The steady version of Eqs. (CD) — (j4j) and boundary conditions @ and (jHJ) 
was discretized by central finite differences. The Euler — Newton and Keller 
[30] arclength continuation algorithms were used [22] to trace out bifurcat- 
ing branches of steady convection using the Harwell MA32 Fortran routine. 
In all these computations, the along-slot period A = 2 was prescribed. The 
grid with 33 nodes in the across-slot direction was used in the computations 
[nx x X(nx + 1) with nx = 33], as in [20.22.24.27J. Temporal behavior of the 
linearized version of Eqs. (CEJ) — (BJ and boundary conditions (jSJ) and (JBJ) was 
also examined, in particular near the onset of oscillatory instability of the con- 
duction state. For this purpose, the implicit method was employed to compute 
time evolution of the respective linear system with time step 5r = 0.05. 

2.2 Linear stability calculations 
2.2.1 Two-dimensional disturbances 

With the state of rest being the background flow for [i = l, the Fourier mode 
of a 2D marginally unstable oscillatory perturbation with angular frequency 
u and wave number k can be written as 



Here [u'(x), t'(x), s'(x)} T is the Fourier-mode part depending on the across- 
slot coordinate alone, the prime near a flow variable denotes such part in 



[u'{x),t'{x),s'{x)] T e iiu)T±ky) +cc. 



(7) 
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the perturbation of the variable. Expression ([Tj) has been introduced into the 
linearized Eqs. (pQ) — (jl]) rewritten for Ra$ = in terms of across-slot velocity 
u, relative temperature t, and solute concentration s. This leads to: 

d 2 d 2 
(^-k 2 )(^- 2 -k 2 - i u J )u = 
dx A dx z 

=F Ra - s) sin 9 ± k 2 (t - s) cos 0] , (8) 
d 2 

(— -k 2 -iuPr)t = u, (9) 
dx 2 

d 2 

(—-k 2 -iuPr)§ = u, (10) 

oar 



where u = —u'Pr, t = t'/Ra, and s = s'/Ra s (Ra s = Ra). 

The variables u, t, and s are subject to the following set of boundary condi- 
tions: 



when the boundaries are stress-free. For no-slip boundaries, the set of bound- 
ary conditions is as follows: 

~ = du =i= d§ =0 (x = ±1/2 ). (12) 
dx dx 



For finding the marginally stable values of Ra and u, it is sufficient to con- 
sider only positive values of u and k in one such traveling wave as ([7]). If 
[u(x),t(x), s(x)] T is a solution of the mode e lt ^ T+kv ^ version of Eqs. (jHJ) — (fTOj) 
at some i?a c and lo c , then the mode e^r-ky) vers i on Q f these equations would 
be satisfied by ±[w(— x), t(— x), s(— x)] T at these same Ra c and uj c . It is thus 
only the mode e l ^ T+hy ^ version of linear stability equations that is hereafter 
implied when such equations are referred to. 

For examination of the linear stability of the conduction state in inviscid fluid, 

d 2 d 
M-r-* - fc2 )« = Ra[ik—(t - s) sin# + k 2 (t - s) cos9] (13) 
dx 2 dx 



was used along with 

(-^-k 2 -iu)t = u (14) 
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and 



"dx 2 



— k — iu)s 



u. 



(15) 



where the Rayleigh numbers are defined as Ra = gaATd 3 / k 2 and Ra s = 
—gP(dS/dx)d 4 /K 2 = Ra. Here u is nondimensionalized with n/d 2 as opposed 
to vjd 2 in Eqs. ([8]) — ([TO]) . Although the definitions of Ra, Ra s , and uj for 
inviscid fluid are different from the respective definitions for viscous fluid, the 
same notations are used for these parameters. It is implied below that the 
definitions of Ra and Ra s as well as the time nondimensionalization scale (in 
uj) correspond to the type of fluid in question. 

The boundary conditions for inviscid fluid are 

u = t = ^ = (a = ±1/2). (16) 



The inviscid problem at 9 = was also examined for the Ra and Ra s being 
independent of each other. That is, equation 

d 2 

iuj(—-k 2 )u = k 2 (Rat-Ra s ~s) (17) 
dx 2 



was considered along with Eqs. (TH|) and (115]) and boundary conditions (fIBl) 
for specified values of Ra s . 

Eqs. ([9]) and (fTTJ]) are first solved for t — s, upon which the resulting general 
solution is introduced into Eq. (jSJ). The latter equation is then solved for u, and 
the general expressions for t and s are thus obtained from Eqs. (jHD and ( TTTJl) . 
For a fixed fc, i?a c and u; c are then found by searching in the Ra — uj domain for 
the smallest Ra at which the complex matrix resulting from the application of 
boundary conditions (TTTT) or (TT2]) to the obtained general solution is singular. 
The same procedure was also applied to the general solution of either Eqs. 
(jT5j) — (IT5]) or Eqs. (JSJ), (TTBT) . and ([17]) along with boundary conditions (]T5"]) . 
NAG Fortran routines were employed for this purpose. 

Once Ra c (ko) and uj c (ko) have been found for a given fe ; the corresponding 
values of these parameters, Ra c (k) and w c {k), at a nearby k = ko + 5k can 
be computed by the Euler — Newton continuation method. This method was 
applied to the solution of equation 

F[Ra(k),w(k),k] = Q, (18) 
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where F[Ra(k),u(k),k] stands for the (complex) determinant of the matrix 
resulting from the application of boundary conditions (TlT!) . (fl2|) . or (fl6l) to the 
general solution of the respective set of differential equations. As mentioned 
above, such determinants were computed with the use of standard Fortran 
routines. Numerical differentiation was thus employed to compute the Jaco- 
bian of {Re[F(Ra, lj, k)], Im[F(Ra, lu, k)]} T (with respect to Ra and lj) and 
dF(Ra,u,k)/dk. 

One also needs to ensure that no disconnected branches of [Ra(k),u(k)] T sat- 
isfying the real and imaginary parts of Eq. (|18p arise for the values of Ra that 
are smaller than those computed by the continuation procedure. For this pur- 
pose, Ra c {k) and u) c {k) were independently obtained for several values of k by 
the search in the Ra — lj domain outlined above. In this case, however, the val- 
ues of k could be relatively scarcely spaced within the considered wave number 
interval. The results of such searches were found to be entirely consistent with 
those obtained by the Euler — Newton continuation. 

For viscous fluid, the linear stability to steady disturbances was also consid- 
ered, as 9 G (tt/2, 7r). For this purpose, lj = was set in Eqs. (JSJ) — (TTDT) and 
the general solution of (resulting) equations 

(^-k 2 ) 2 u = -Ra[ik-^-(t-§) sm6 + k 2 (t-S)cose], (19) 
ax 2 ax 



and 



"dx 2 



'dx 2 



-k 2 )t = u, (20) 



k 2 )s = u (21) 



was obtained analytically. Either boundary conditions (fTTj) or 0121) were then 
applied to this general solution and the smallest Ra, Ra c (k), at which the 
resulting matrix becomes singular were searched for at different k. 



2.2.2 Three-dimensional disturbances 

For 9 ^ 0, 7r, a three-dimensional (3D) structure of perturbation can generally 
become relevant in the linear stability analysis. This study is mainly focused on 
2D disturbances, for which understanding of the basic physical effects of differ- 
ent boundary conditions can be most conveniently obtained. Unless explicitly 
stated otherwise, therefore, it is the 2D linear instability that is implied below. 
However, 3D disturbances were also considered: in some cases, their analysis 
was found to be useful for understanding the 2D results. 
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With the same nondimensionalization as in u, t, and s, the Fourier mode of a 
3D disturbance is described as 

[u(x),t(x), S ( x )f e ^T+k y y+k z z) + (22) 



where k y and fc z are the y and z components, respectively, of the full wave 
number k = (k 2 + k 2 ) 1 ^ 2 . (The 2-axis is orthogonal to the x — y plane in Fig. [T] 
and is directed towards the reader.) Expression ff22l) was then introduced into 
the appropriately nondimensionalized, linearized 3D governing equations for 
viscous and inviscid fluid. 

With the definition of wave number k just generalized for 3D disturbances, the 
marginal linear stability in viscous fluid is described by the equations formally 
identical to Eqs. ([9]) and (TTUT) along with 

d 2 d 2 „ d 

^dx 2 ~ ^^dx 2 ~ ^ ~ %UJ ^ = ~ Ra][lkv di^~ ^ Sin ^ + 

k 2 (t- 5) cos 9]. (23) 



Likewise, the 3D marginal stability in inviscid fluid is described by the equa- 
tions formally identical to Eqs. ( fT4l) and ( fT5l) along with 

d 2 d 
icu(— - k 2 )u = Ra\ik v — (t - s) sin6> + k 2 (t - s) cos 01. (24) 
dx 2 y dx 

For fixed values of k z , Ra c (k y ) and u> c (k y ) were obtained with the same pro- 
cedure as described above for 2D disturbances. 



3 Oscillatory convection 



3.1 General 



For interpretation of the findings on oscillatory instability, only a reflectionally 
symmetric oscillatory perturbation is considered below. In a horizontal slot, 
such perturbations are often referred to as standing waves. They are then char- 
acterized by convection cells changing their sense of rotation periodically in 
time. As suggested by [31], both reflectionally and translationally symmetric 
oscillatory perturbations are expected to arise from the Hopf bifurcation in 
a system whose background state possesses both these symmetries. Transla- 
tionally symmetric oscillatory perturbations are often referred to as traveling 
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waves. However, it is sufficient to understand the physical mechanism of such 
oscillatory instability in terms of only one of the types of perturbation just 
mentioned. Amplitude growth of the oscillatory perturbation of the other type 
can then be viewed merely as a mathematical consequence of the respective 
results in [3T] . 

3.2 Inviscid fluid 

3.2.1 Effect of the across-slot gravity. 

Let us first consider the problem in a horizontal slot addressed in [19J (6 = 0), 
whose linear stability is described by Eqs. (JHJ), ffl5]) . and (ITT)) and boundary 
conditions (fTBj) . The marginal-stability curves, Ra c (k) and oo c (k), are illus- 
trated in Fig. [2j They exhibit two basic features. The first is that the most 
unstable wave number is zero for any Ra s . This wave number is also char- 
acterized by zero oscillation frequency. The second feature is the existence of 
a minimal horizontal wavelength, decreasing with the increase of Ra s , below 
which the instability does not arise. 

Reinterpreting [19], rotation of a small-amplitude perturbation cell generates 
a potential energy of component perturbation stratifications in the end of a 
rotation cycle. Due to differential gradient diffusion, this potential energy is 
utilized by the perturbation cell in the beginning of the cycle of cell rotation 
in the opposite sense. The maximal amount of such energy depends on the 
time available for a fluid element of the cell to change its vertical coordinate. 
This time is specified by the horizontal scale of the instability. For this reason, 
Ra c (k) decreases with the increase of the horizontal wavelength and assumes 
its minimal value for a given Ra s as k — ► (i.e., as wavelength A — > oo). 

Below a certain critical wavelength, the potential energy utilized by a pertur- 
bation cell becomes insufficient for bringing about growth of the perturbation 
amplitude. The increase of Ra s and associated growth of Ra c enhance the dis- 
parity between the respective gradients in the perturbed state. This makes the 
energy transfer to the perturbation cells more intensive. The minimal unstable 
wavelength thus decreases as Ra s grows (Fig. [2]). 

Efficiency of the utilization of the potential energy by a perturbation cell de- 
pends on the frequency with which the marginally unstable cells change their 
sense of rotation. This frequency thus has to be such that the time for ver- 
tical diffusion naturally specified by the instability wavelength be resonantly 
matched. Since such time grows with the wavelength increase, u c (k) — > as 
k — > and u c (k) grows with the increase of k from 0. As the wave number 
is further increased, the wavelength time for the manifestation of differential 
diffusion eventually becomes insufficient for the cell oscillation amplitude to 
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grow. As a consequence, u c (k) begins to decrease when certain values of k are 
exceeded (Fig. [2]), to afford more time for the diffusion. 

The decrease of u c (k), however, leads to an inconsistency between the diffu- 
sion time afforded by the perturbation frequency and that specified by the 
instability wavelength. This results in the efficiency of utilization of the com- 
ponent potential energy by a perturbation cell being reduced. The additional 
energy that can actually be utilized due to such frequency decrease is thus ex- 
pected to be limited. As the wave number exceeds a critical value, therefore, 
the instability fails to develop. 

The fact that the most unstable wave number is zero makes the determination 
of exact values of Ra c (0) and group velocity a>£(0) = du) c (0)/dk relevant. Also 
applicable to Eqs. (fT51) — (fT5|) for 6 > 0, the long-wavelength expansion used 
for this purpose is as follows: 

Ra c (k) = Rao + k 2 Ra 2 + 

uj c (k) = k(u)o + k 2 uj 2 + ...), 

u(k) = k(uo + ku\ + k 2 u 2 + •••)> 

t(k) =t + kti + k 2 t 2 + 

s(k) = so + Hi + k 2 § 2 + .... (25) 

Upon introduction of (1251) into Eqs. (tHJ) . (TT5l) . and ifTTj) . and use of boundary 
conditions f|T6l) . to, So, and uo are first obtained. The k 1 order of Eq. ffT5l) and 
the flux-free boundary conditions for s\ from ffT6|) then yield 



Having derived the expressions for s\ and i\ and on introduction of these into 
the k 1 order of Eq. (fl71) . one obtains the expression for u\. With such Ui, the 
k 2 order of Eq. ffT5|) and the flux-free boundary conditions for s 2 from (JT6]) 
thus yield 



The numerical data underlying the marginal-stability curves in Fig. [2] were 
found to accurately coincide with Eqs. (1261) and ( 1271) . For the case of two 
compensating gradients this work is aimed at, Ra c (0) = Ra s , 




(26) 



Ra c (0) = Ra = (2Ra s + 5040) /51. 



(27) 




(28) 
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3.2.2 Effect of the along-slot gravity. 



For the horizontal and inclined slots with two compensating across-slot gradi- 
ents, Ra c {k) and LU c (k), obtained from Eqs. ([TBI — (fT5|) and boundary condi- 
tions (|T6|) . are shown in Fig. [31 (In this Sec. 13.2.21 k z = and thus also k = k y 
are implied.) As in the horizontal slot [Figs. [2] and 0(a)], = remains the 
most unstable wave number up to the orientation of the slot being nearly ver- 
tical [Figs. MJo) — E^ e )]- With increasing 9, however, Ra c {k) and ui c {k) mostly 
decrease and qualitative changes arise in the shape of these curves [Figs. [31(e) 
and[3jf)]. This suggests the emergence of an additional factor destabilizing 
the system. 

As suggested by [2TJf2"Tj . if 9 = a viscous fluid in the region of Ra < 
is expected to develop a steady instability whose mechanism is conceptually 
analogous to the finger instability in the classical double-diffusion |16j . Such 
mechanism is the only cause of linear steady instability for < 9 < it/ 2 
as well (see Sec. 14.1.11 below). This steady instability would thus have finite 
negative values of the marginally unstable Rayleigh number for viscous fluid 
[2"Tj (see also Sec. 14.1.21) . Eqs. fiBfy — (J2Tj) then suggest that in the inviscid fluid 
(Pr — > 0), Ra c (k) = for the respective wave number interval. Ra < would 
thus be the unstable region, with the instability mechanism being similar to 
that discussed in [2"0f2~T] . Ra > is therefore the region of stability to the 
effect of the across-slot gravity component on steady disturbances. 

As the along-slot gravity component arises with 9 growing from 0, however, 
another mechanism of differential gradient diffusion becomes increasingly rele- 
vant. This mechanism favors a monotonic growth of disturbances. For viscous 
fluid, it has been discussed in [221124] . In the absence of dissipation (through- 
out this work, this term implies only viscous dissipation), such mechanism is 
reasonably expected to be effective for small-amplitude disturbances at any 
Ra 0. (The results in Fig. [3] discussed just below confirm this assumption.) 
So long as 9 < 7r/2, however, it should not give rise to steady linear instability 
for Ra > 0. The differential diffusion due to the across-slot gravity component 
opposes amplitude growth of the perturbation cells whose sense of rotation 
does not change (Sec. 13.2. ip . 

For < 9 < 7r/2, part of the rotation energy of a perturbation cell thus comes 
from the energy directly contributed via the along-slot gravity component in 
the current cycle of rotation [2~2~f2~4] . In the end of a cell rotation cycle, the 
whole rotation energy is transformed into the potential energy of the pertur- 
bation stratification due to the across-slot gravity component. This potential 
energy is released in the next rotation cycle. Depending on the relative roles 
of the across-slot and along-slot gravity components, u c (k) thus decrease with 
increasing 9. 
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For Ra cos 9 being between and the Ra c (k) from Fig. [3ja), the growth of 
oscillatory perturbations is largely due to the monotonic contribution of the 
along-slot gravity component, effective at any Ra 7^ 0. The relative role of the 
across-slot gravity component in such growth on a given wavelength is specified 
by the respective effectiveness of transformation of the energy of cell rotation 
into the potential energy of the across-slot perturbation stratification. (Such 
process maintains the oscillatory nature of the instability.) This effectiveness 
can be judged based on the behavior of Ra c {k) in Fig. [3j(a). 

Compared to the vicinity of the infinite wavelength, finite wavelengths are less 
effective in transforming the energy of cell rotation into the above potential 
energy. This makes the relative role of the along-slot gravity component on 
such finite scales more pronounced than for the longer wavelengths. [As 9 grows 
from in Figs. [3(a) — [3J(d), in particular, Ra c {k) decrease for k >~ 1.1 and, 
due to the decrease of the across-slot gravity component, increase slightly for 
k <~ 1.] Near 9 = it/2, a set of finite wavelengths thus becomes more unstable 
than the vicinity of the infinite wavelength [Fig. [3J(f)] . 

With the across-slot gravity component vanishing as 9 — > n/2 [Figs. [3(e) and 
[3(f)], however, both Ra c (k) and uj c (k) decrease for all k. For 9 > ir/2, therefore, 
the slot with inviscid fluid is expected to be unstable to a continuum of the 
wavelengths of steady disturbances for any Ra > 0, and it ought to be unstable 
for any Ra 7^ at 9 = tt/2. As mentioned above, this would be consistent with 
the interpretations in [221T21j if the absence of dissipation is allowed for. The 
oscillatory and steady instabilities due to the respective effects of differential 
gradient diffusion in inviscid fluid are thus continuously transformed into each 
other around 9 = it/ 2. 

Upon introduction of fT25]) into Eqs. ([TBI — (TT5"j) . the k 1 order of these equations 
and boundary conditions ([16]) yields 



Ra 2 sin 2 9 + 720Ra cos 9 - 8640cj 2 = 0, (29) 



which is consistent with ( |26|) for 9 = 0. The values of Ra c (0) = Rao and 
u>l(0) = ujq estimated from the numerical data underlying the (2D) marginal- 
stability curves in Fig. [3jwere found to be consistent with Eq. (!29j) to a fairly 
high accuracy. 
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3.3 Viscous fluid 



3.3.1 Effect of the across-slot gravity. 

The curves of 2D marginal stability for viscous fluid under stress-free and no- 
slip boundary conditions are illustrated in Figs. H] and [5j They were obtained 
from Eqs. (jHJ) — (fTOj) and boundary conditions ffTTl) and ffl2l . respectively (In 
this Sec. 13.3.11 k z = and thus also k = k y are implied.) The most prominent 
feature of the 2D curves in Figs. H] and \5[ compared to the above results 
for inviscid fluid, is that k = is not the most unstable wave number and 
u> c (k) > as k — > 0. 

In the context of the effect of the across-slot gravity component, the obtained 
stability of the infinite wavelength may at first sight appear to be counterin- 
tuitive. The overall dissipation of perturbation motion is minimized when the 
wavelength becomes infinite. Since the effect of differential gradient diffusion 
is maximized for A — > oo, one may think that the infinite wavelength has to 
remain most unstable in viscous fluid as well. In particular, k = is most un- 
stable to steady disturbances in the viscous fluid problem identified by 9 = ir 
and Ra > [23ME2]. 

The relative stabilization of the long wavelengths [seen in Figs. H](a) and[5](a), 
for example] with respect to the action of the across-slot gravity component is 
associated with a special feature of the effect of dissipation. This feature arises 
when the instability is oscillatory. Before describing such feature, however, it 
is worth illustrating details of the oscillatory instability mechanism in the 
horizontal slot. 

Fig. |6J(a) represents the stage when the previously accumulated potential en- 
ergy of perturbation stratification is released after the sense of cell rotation has 
changed. At this stage, cell motion is favored and opposed by the horizontal 
solute and temperature perturbation gradients, respectively. Due to differen- 
tial (across-slot) gradient diffusion, the temperature opposition to convective 
motion in Fig. [6]^a) is overcome. Such motion thus intensifies [see ip in Figs. 
E(a) andEKb)]. 

At the current stage of the oscillation period, the release of the accumulated 
potential energy is effected via a mechanism similar to that discussed in [20||2~T] . 
where an inversely stratified problem was considered. This is particularly ap- 
parent as one compares the present Fig. El^a) with Fig. 2(a) in [21]. In the 
present problem, however, the intensification of convective motion [tp in Figs. 
[6ja) and[6j(b)] is inevitably accompanied by changing distributions of the com- 
ponent perturbation isolines [s and t in Fig. [61(b)] . With the lower background 
values of the components at the upper wall, the clockwise (counterclockwise) 
rotation gives rise to the relatively smaller respective perturbations in the 
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right (left) part of a rotating cell. The component perturbation isolines are 
thus gradually reset so that cell motion be favored and opposed by the tem- 
perature and solute perturbation stratifications, respectively [Fig. 0(c)]. 

Due to the interaction of a rotating perturbation cell with flux boundary 
conditions, the solute perturbation isolines become more and more vertically 
oriented, in contrast to the perturbation isotherms [s and t in Figs. [6]^d) and 
E(e)]. The resulting differential gradient diffusion acts against the sense of 
rotation of a convective cell. This damps cellular motion [Figs. EIc) — El(e)]. 
The energy of cell rotation is thus transformed into the potential energy of 
the component perturbation stratifications. Such potential energy is released 
[Figs. [6](f) andH^g)] after the sense of cell rotation has changed. 

The events in Figs. El^f) — Mj) are thus qualitatively identical (up to the inverse 
perturbation sign) to those in Figs. 0(a) — El^e). The streamline pattern in Fig. 

is different from that in Fig. [6](e). This is associated with the stages of the 
change in the sense of cell rotation being slightly different in these figures. Fig. 

thus shows that the change of rotation sense is initiated via the formation 
of small cells of the opposite sense of rotation near the boundaries, where 
differential gradient diffusion is most effective. 

The oscillatory instability thus arises when the potential energy generated 
during the rotation of a perturbation cell in one sense is sufficient to increase 
the cell amplitude while the energy is released. The release of potential energy 
is effected via differential gradient diffusion during the cell rotation in the 
opposite sense. Part of such energy, however, is inevitably spent on dissipation 
of along-slot cell motion. This part is the greater the larger the wavelength 
is. Efficiency of the (oscillatory instability) feedback is thus reduced by the 
dissipation. [Such feedback forms between the perturbation gradient disparity 
in one cycle of cell rotation and the amplitude of rotation (giving rise to such 
a disparity) in the next cycle.] The feedback efficiency, therefore, also tends 
to zero when the wavelength becomes infinite. 

Although increase of the instability wavelength decreases the overall dissipa- 
tion, it enhances the role of the unchanged part of the dissipation (i.e., of the 
dissipation of along-slot cell motion) in the oscillatory instability feedback. 
For sufficiently long wavelengths, the latter factor dominates the overall effect 
of dissipation. To avoid the zero feedback efficiency when oj c — > 0, u) c {k) thus 
remains greater than zero as k — > 0. In the long- wavelength limit, therefore, 
an increasing portion of fluid particles undergoing a cycle of differential gradi- 
ent diffusion fail to acquire a sufficient potential energy for their perturbation 
amplitude to grow in the next cycle. In the context of the effect of the across- 
slot gravity component, this explains the relative stabilization of long (2D) 
wavelengths in Figs. H](a) — @Jd) and [5]^a) — [5](d) compared to the respective 
inviscid-fluid problems [Figs. [3(a) — [21(d)] . 
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3.3.2 3D effects. 



The above interpretation of the effect of viscous forces on the oscillatory insta- 
bility suggests that three-dimensionality of the disturbances might be relevant 
in an inclined slot with viscous fluid. Let the 2D perturbation wavelength, X y , 
be large enough for the effect of dissipation on efficiency of the instability 
feedback to dominate the overall effect of dissipation. Such efficiency increases 
as the orientation of the axis of rotation of a convective cell changes and the 
total wavelength decreases with the emergence of a z component of the wave 
number: A = 2ir/(k 2 + k 2 ) 1 / 2 < \ y = 2ix/k y . The critical Rayleigh number 
could thus decrease compared to the 2D problem. 

The horizontal-slot linear stability problem depends only on the wave number 
modulus. This implies that 

Ra c (k y , k z ) = Ra c [{k 2 + k 2 z ) 1/2 , 0]. (30) 



Let k y ( > 0) be fixed, with respect to the selected coordinate system, in the 
interval where Ra c (k y )\k z =o decreases with increasing k y . Then Ra c {k y )\k z >o < 
Ra c {k y )\k z =o for some k z > 0. That Ra c (k y )\ kz>0 > Ra c (k y )\ kz=0 for all positive 
k y and k z in inviscid fluid also follows from the Ra c {k y ) in Fig. [3](a) only 
growing with increasing k y . 

When 9 > 0, the wave number orientation becomes physically meaningful. 
For small 9 > 0, however, Ra c (k y )\k z >o are still largely determined by the 
linear stability equations for 9 = 0. These describe the zero-order perturbation 
expansion in 9 of such equations for 9 > 0. When 9 is small, therefore, most 
unstable disturbances with small k y are still expected to be of 3D and 2D 
nature in viscous and inviscid fluid, respectively. 

Such a three-dimensionality of instability hinges on the existence of an inter- 
val of growing k with decreasing Ra c (k) in a parametrically close system of 
marginal-stability equations for which Eq. fl30l holds. It could thus also apply 
to different problems, where such formal conditions for the respective stability 
parameter are met due to other physical effects. 

For a finite 9 > 0, the role of 3D disturbances is no longer formally predictable. 
However, the above physical interpretation of the effect of dissipation on the 
instability feedback still suggests the 3D nature of the most unstable distur- 
bances with relatively small k y . This effect is also due to the across-slot gravity 
alone. It therefore has to vanish as the slot orientation approaches vertical. 

Examination of 3D oscillatory disturbances in the inviscid fluid [Figs. [3](b)— 
[3]^e)] shows that the disturbances with k z up to 3 are more stable than the 
respective 2D perturbations [k z = 0). Meaning that small- amplitude 3D dis- 
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turbances of any physical nature are not expected to arise in the range of 
parameters considered in Figs. [3(b) — [3(e), this result is consistent with the 
interpretation just mentioned. 

As also anticipated from the above interpretation, 3D oscillatory disturbances 
are most unstable in viscous fluid for an interval of k y adjoining k y = 
[Figs. [4(b)— i(e) and 0(b)— [5(e)] . In addition, the data in Figs. Kb)— He) 
and [5(b) — [5(e) indicate that 3D disturbances are the first to arise only when 
the across-slot gravity is present. [2D disturbances {k z = 0) are most unstable 
for all k y presented in Figs. H(e) and 0(e).] Consistent with the above inter- 
pretation of the effect of the across-slot gravity component, Figs. H(e) and 
E(e) thus also suggest that the along-slot component does not introduce an 
additional 3D oscillatory instability. 



3.3.3 Effect of the along-slot gravity. 

3.3.3.1 General. As the results in Figs. H(e) and [5(e), apart from the 
respective findings in [21], suggest, 2D oscillatory instability also arises due to 
the effect of the along-slot gravity component alone. The mechanism of such 
instability can be identified by careful examination of Fig. In this figure, 
temporal behavior of an unstable perturbation mode possessing one of the 
reflectional symmetries of the linearized problem is illustrated. 

The perturbation flow structures illustrated in Fig. [7] suggest that their dynam- 
ics is underlain by two counter-propagating patterns. Either such pattern is 
dominant near one of the slot sidewalls. The pattern propagating in (against) 
the gravity direction, i.e. downwards (upwards), is dominant near the left 
(right) sidewall in Fig. This could be viewed as manifestation of such a per- 
turbation being a superposition of two counter-propagating traveling waves. 
In particular, let g(r, x, y) stand for the reflectionally symmetric perturba- 
tion of if>, t, and s, nondimensionalized say consistently with Eqs. ([T]) — (jlj). If 
Ra c {k) < Ra c {nk) (n = 2,3,...), in the large-r marginally unstable state of 
such a wave number k, 

q{r,x,y) = Re[a\q' {x)e^ cT+ky) + a q 2 q'{-x)e i{ ^ T - ky) ] 
= ReiaWixy^+^il + 4^f^e- 2 ^]} 

= Re{a q 2 q'(-x)e l ^ T - k ^[l + ^-^L e 2iky }}. (31) 

Here coefficients a\ and a\ specify the symmetry of the respective variable, and 
q'(x) is the part of the Fourier mode of the variable depending on x alone. It 
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is seen from Eq. (I3T]) that when <C x)| for x > 0, the slot half near 

the left (right) sidewall is dominated by the downwards-(upwards-)propagating 
traveling wave. Dynamics of such a reflectionally symmetric mode restricted 
to a slot half is thus largely due to the dominant part of the respective prop- 
agating pattern. 

Such dominant parts of the downwards- and upwards-propagating stream- 
line perturbation patterns are distinguishable as the along-slot sequences of 
counter-rotating cells in the left and right halves, respectively, of Figs. [7(d) and 
[3^i)- Near the left (right) sidewall, the regions of the clockwise- rotating cells 
from such along-slot cell sequence are seen from Fig.[7]to propagate downwards 
(upwards) concurrently with mostly negative (positive) perturbation values of 
both components. The regions of such counterclockwise-rotating cells are also 
seen to propagate downwards (upwards) concurrently with mostly positive 
(negative) component perturbation values. Why such component perturba- 
tion distributions take place is analyzed below. 



3.3.3.2 Component perturbation distributions. Let us consider the 
dominant part of either traveling streamline perturbation pattern [Figs. [3(d) 
and [3(i)] in the respective frame of reference moving with it. In the refer- 
ence frame moving with the downwards-(upwards-)propagating pattern, the 
laboratory-frame velocities (Fig. [7j) oriented downwards (upwards) decrease. 
Such velocities oriented upwards (downwards) thus increase. For an adequate 
speed of propagation, the component perturbations are therefore transported 
vertically only against the direction of propagation of the respective pattern. 

The distributions of temperature and solute perturbations in Figs. [T^d) and 
[3^i), in particular, are in part the result of the vertical transport just high- 
lighted. They are also due to the advection associated with across-slot laboratory- 
frame velocities. The orientation of such an advection is the result of addition 
of the along-slot component of motion against the respective direction of prop- 
agation to these laboratory-frame velocities. 

Allowing for the along-slot periodicity, the regions near the upper and lower 
domain boundaries in Fig. 0^d) are in the wake of the downwards-propagating 
clockwise-rotating cells. Such a cell is seen in the left central part of the flow 
domain in this figure. These regions are also in the wake of the upwards- 
propagating counterclockwise-rotating cells represented by such cell in the 
right central part of the flow domain [Fig. [3(d)] . The component perturba- 
tions are thus transported into the upper- and lower-boundary regions by 
the velocities whose (laboratory-frame) orientation in such convective cells is 
against the respective direction of propagation. 

The component perturbations in the upper- and lower-boundary regions in 
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Fig. [7(d) are also affected via the advection caused by (laboratory-frame) 
across-slot velocities. On addition of the along-slot component of motion against 
the respective direction of propagation, such across-slot velocities in the left 
(right) half of the slot are relevant only where they are oriented in Fig. [3(d) 
from (towards) the respective sidewall. 

In all transport processes just outlined, the temperature and solute perturba- 
tions are transferred into the upper- and lower-boundary regions [Fig. [7(d)] 
only from the near- left-wall and middle areas of the slot. In such areas, the 
background values of the components are in the upper half of their interval. 

The central region in Fig. ETJd), occupied by the full counter-rotating cells, is 
also in the wake of two counter-propagating convective cells. One of these is 
the upwards-propagating clockwise-rotating cell partly seen in the right upper 
corner of the domain in Fig. [7Jd). The other is the downwards-propagating 
counterclockwise-rotating cell partly seen in the left lower corner of the do- 
main. The component perturbations are thus transported into the central 
region of the domain in Fig. [3(d) by the velocities whose (laboratory-frame) 
orientation in such convective cells is against the respective direction of prop- 
agation. Advection due to the (laboratory- frame) across-slot velocities just 
below (above) the domain central region in the left (right) slot half in Fig. 
[TJd) also affects the component perturbations in this region. 

In all transport processes just outlined, the temperature and solute perturba- 
tions are transferred into the central region of the domain in Fig. [3(d) only 
from the near-right-wall and middle areas of the slot. In such areas, the back- 
ground values of the components are in the lower half of their interval. 

The central region of the domain in Fig. [7J(d) is thus characterized by the 
relatively smaller values of the component perturbations while the larger val- 
ues are in the upper and lower regions of the domain. The dominant parts 
of counter-propagating streamline perturbation patterns distinguished in Fig. 
[TJi) are about half of the wavelength ahead of their locations in Fig. [7Jd) in 
the respective directions. The component perturbation distributions in Fig. 
[TJi) are consistent with the above understanding [described for Fig. [3(d)] of 
the transport processes given rise to by the cells in this figure. 

The component perturbations in Fig. [7J(i) can also be viewed as resulting 
from their dominant parts being about half of the wavelength ahead of their 
locations in Fig. [TJd) in the respective directions. These dominant parts are 
however indistinguishable in Figs. [3(d) and[7Ji), where they merge into the 
respective slot-wide patterns of the component perturbations. They are only 
approximately discernible in Figs. [TJa) and[7Je). 

As co-rotating cells in the dominant parts of two counter-propagating stream- 
line perturbation patterns [Figs. [TJd) and [7J(i)] have close vertical locations, 
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they merge into a single slot- wide cell. This is seen in Figs. [7(a) — E[c) j7(e)J3(f)— 
[TJh), and[3(j)- The temperature and solute perturbations in these figures result 
from superposition of their respective counter-propagating perturbation pat- 
terns [whose dominant parts are roughly discernible in Figs. 0(a) and[3(e)]. The 
latter accompany the counter-propagating streamline patterns [whose domi- 
nant parts are discerned in Figs. [3(d) and[?Ji)]- 

The component perturbation distributions in Figs. [3(a) — [7(c) j7Je)J3(f) — E(h), 
and[3j') are thus also consistent with the transport processes given rise to by 
the respective patterns of convective cells in these figures. Such processes have 
to be interpreted as in the above analysis of Figs. [3(d) and[7Ji)- 



3.3.3.3 Feedback. In the pattern dominant part propagating downwards 
(upwards) near the left (right) sidewall, the clockwise-rotating cells thus steadily 
coincide with the regions of mostly negative (positive) perturbation values of 
the components. Such counterclockwise-rotating cells then steadily coincide 
with the regions of mostly positive (negative) component perturbation values 
[Figs. [3(d) and[7Ji)]- In the vicinity of a sidewall, the perturbation isotherms 
are however prone to be parallel to the boundary, where the temperature per- 
turbation vanishes. In contrast to this, the solute perturbation isolines tend to 
be orthogonal to the sidewalks, because of the flux-free boundary conditions. 

Such disparate isoline behavior gives rise to differential gradient diffusion. As a 
consequence, a near-slot-middle streamline point in the downwards- (upwards- 
)propagating clockwise-rotating cells [Figs. [3(d) and[7Ji)] is typically heavier 
(lighter) than the respective streamline point near the left (right) sidewall at 
the same vertical location. Likewise, a near-slot-middle streamline point in the 
downwards- (upwards-)propagating counterclockwise-rotating cells is typically 
lighter (heavier) than the respective streamline point near the left (right) 
sidewall at the same vertical location. [In particular, see Figs. [8(d) and[8](i).] 

The flow perturbation structures illustrated in Fig. [3 are expected to result 
from superposition of the counter-propagating patterns whose dominant parts 
have just been described. They are all thus characterized by a horizontal den- 
sity difference between two streamline points that is consistent with the sense 
of rotation of the respective convective cell. This can be inferred from Fig. [3 
and is explicitly illustrated in Fig. [SJ 

In a quasisteady sense, therefore, the energy that drives small-amplitude os- 
cillatory convection in Fig. [3 has essentially the same nature as that driving 
small- and finite-amplitude steady convection in [22.24J . In contrast to small- 
amplitude steady convection, such oscillatory convection arises for fi = 1 be- 
cause the oscillatory perturbation gives rise to a more favorable ratio between 
the maximal absolute values of the relative perturbations of the components. 
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This is discussed in Sec. 14.1.11 below, in the context of steady convection. 



The streamline disturbance resulting from superposition of certain counter- 
propagating convective flow patterns thus generates such distribution of the 
component perturbations as favors its convective motion. For amplitude growth 
of such a disturbance to be maintained, however, the counter-propagation of 
the convective patterns also has to be sustained. This underlies the oscillatory 
nature of the instability. 

It was assumed above that, in the marginally unstable state, the speed of 
streamline-pattern propagation adequately matches the laboratory-frame con- 
vective velocities oriented in the respective direction of propagation. The cor- 
responding streamline points thus have their component perturbation values 
practically steady. 

In the counter-propagating sequences of convective cells [Figs. [3(d) and [31)], 
practically steady values of the component perturbations arise in the near- 
wall regions of the counterclockwise-rotating cells. The values of component 
perturbations in these cell sequences also have to be practically steady in 
the middle slot regions of the clockwise-rotating cells. For all these regions, 
the laboratory-frame velocities are oriented in the direction of propagation 
of the respective streamline pattern. In the middle slot regions, the effects 
of temperature and solute perturbations on the density perturbation largely 
offset each other. Due to differential gradient diffusion, however, the density 
perturbation in the near-wall regions is specified primarily by the respective 
values of solute perturbation alone. 

The solute (density) perturbation values in the near-wall regions of counterclockwise- 
rotating cells in Figs. [7(d) and [7J(i) [Figs. [8(d) and Mj)] are expected to be 
steadily positive and negative in the left and right parts of the slot, respec- 
tively. The time instances when such counter-propagating counterclockwise- 
rotating cells merge into a single counterclockwise-rotating cell are illustrated 
in Figs. [7(a) — [7J(c)j7(e)|7J(f) — [7(h), and UJj). The values of solute (density) 
perturbation in such a slot-wide cell near the left sidewall are also always 
larger than those near the right sidewall (Fig. [8]). 

It is the steady horizontal density difference arising between two halves of the 
slot that sustains the counter-propagation underlying dynamics of the reflec- 
tionally symmetric disturbance. The counter-propagating convective stream- 
line patterns thus give rise to such feedback distribution of the component 
perturbations as, due to differential gradient diffusion, enhances their convec- 
tive motion and sustains their counter-propagation. 
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3.3.3.4 Marginal-stability boundaries. Differential gradient diffusion 
is most effective when a convective cell particle has the maximal time available 
at a particular wavelength to unequally change its component values while its 
across-slot location is changed. This time is the larger the longer the wave- 
length is. In order that the oscillatory instability mechanism just described be 
effective, such diffusion time has to be adequately matched by the time during 
which the wavelength is passed by either traveling perturbation pattern. The 
latter time specifies the marginal instability frequency. Such frequency in Figs. 
@Je) and [Sf e) thus decreases as the wavelength increases. (In this Sec. 13.3.3.41 
k z = and thus also k = k y are implied.) 

Dissipation of the along-slot component of perturbation motion is part of the 
feedback for the steady instability caused by different boundary conditions 
in a vertical-slot geometry (see Sec. 14.1.11 below). The efficiency of such a 
feedback thus has to decrease as the wavelength increases, as discussed for the 
horizontal-slot (9 = 0) oscillatory instability feedback (Sec. I3.3TT]) . 

In view of the vertical-slot oscillatory instability feedback being of a quasis- 
teady nature, its efficiency is also progressively reduced by the dissipation 
of along-slot perturbation motion when the wavelength increases. As in the 
horizontal-slot oscillatory instability, the feedback efficiency would also de- 
crease to zero if the time available for across-slot diffusion as A — > oo {k — > 0) 
were fully utilized. In the long-wavelength limit, therefore, u c (k) in Figs. @Je) 
and [3^e) remains a finite nonzero value. However, differential gradient diffu- 
sion, being a basic ingredient of the vertical-slot oscillatory instability mech- 
anism, would then have to progressively decrease its effectiveness as k ap- 
proaches 0. This explains why Ra c (k) increase as k decreases near k = in 
Figs. H](e) and [5](e) and the (2D) marginal-stability curves in these figures 
[as well as those in Figs. @](b) — S](d) andE^b) — E3(d)] remain qualitatively the 
same as in Figs. S](a) and 0(a). 

As 9 increases above it/2, the emergence of the across-slot gravity component 
makes the associated mechanism of linear steady instability [201121] potentially 
relevant. For a particular wavelength, the diffusion time thus has to optimally 
match a combination of the steady manifestation of the effect of the across-slot 
gravity component and the oscillatory one of the along-slot component. The 
marginally unstable frequency, uj c (k), is thus generally expected to decrease as 
9 exceeds n/2. As the wavelength increases, however, such frequency specified 
by the effect of the along-slot gravity component alone progressively exceeds 
that corresponding to the maximal time available at the wavelength for the 
respective manifestation of differential gradient diffusion. This allows to avoid 
the overly large reduction of the feedback efficiency. 

So long as the effect of the along-slot gravity component is dominant, the 
dissipation of along-slot perturbation motion remains of primary importance. 
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Intermediate and long wavelengths could not thus significantly decrease their 
marginally unstable frequency. For this reason, u c (k) decreases noticeably only 
in the region of relatively large k in Figs. H(f) and [5(f), whereas it remains 
practically unchanged in the vicinity of /c = 0. The intermediate and long 
wavelengths thus fail to substantially absorb the steady manifestation of the 
effect of the across-slot gravity component. It is only for the relatively large 
wave numbers, therefore, that Ra c (k) decrease with increasing 9. In contrast 
to this, the decreasing effect of the along-slot gravity component leads even 
to slight increases of Ra c (k) in the interval of relatively small wave numbers 
in Figs. Bp) and [5(f). 

With 9 increasing above n/2, the decrease of the along-slot gravity component 
still enhances room for the growing across-slot component. The lower bound of 
the region of k with decreasing Ra c (k) is thus found to decrease [Figs. [4(f) and 
[5(f)]. Around 9 w 1.57r/2, however, Ra c (k) are still increasing (with increasing 
9) when k <~ 2.5 for both stress-free and no-slip boundary conditions. In 
particular, the development of a new minimum in the shape of Ra c (k) [Figs. 
H]T ) and [5(f)] is a manifestation of the opposing trends exhibited by Ra c (k) 
within the different parts of the interval of k. 

For 9 ~ 1.557r/2, u c (k) in the regions of the smallest and largest values of k in 
Figs. H]T ) and [5(f) could not be continued into each other. Slope \duj c (k)/dk\ 
also seemed to exhibit infinite increases around an interval of middle values 
of k. No marginally unstable curve with u c {k) > was then found within 
this interval of middle k. These findings suggest that at 9 « (1.5 — 1.55)7r/2 
the oscillatory marginal-stability boundaries in Figs. H(f) and [5(f) begin to 
transform into the respective curves characterized by a zero real eigenvalue. 



4 Steady convection 

4-1 Linear steady instability 

4-1.1 Effect of the along-slot gravity. 

The vertical slot with inviscid fluid was found in Sec. l3.2.2l above to be unstable 
to steady infinitesimal disturbances for a continuum of wave numbers at any 
Ra 7^ 0. In viscous fluid, however, dissipation of the along-slot component of 
perturbation-cell motion is part of the feedback for steady instability, if any, in 
the vertical-slot geometry. In particular, the component gradient disparity and 
resulting differential gradient diffusion are caused by across-slot perturbation 
motion. Such across-slot motion is however maintained only due to along- 
slot motion in the convective cell. Driven by the horizontal density differences 



25 



between respective streamline points, such along-slot cell motion arises from 
differential gradient diffusion [22.24J, the process responsible for the feedback. 



For viscous fluid, the linear steady instability of all wavelengths is thus gener- 
ally less likely in a vertical-slot geometry than in the horizontal slot where the 
component stratification favors steady instability [2"0|2~T] . In addition, if such 
a vertical-slot instability does arise (as in [22], for example, where the back- 
ground horizontal gradients are not exactly compensating), the long wave- 
lengths are not expected to be most unstable. This is associated with the 
feedback efficiency tending to zero as A — ► oo. As seen from Figs. M arid [TO} 
however, the long wavelengths are most unstable for all illustrated 9. Their 
instability is also most persistent as the interval of unstable wave numbers 
shrinks and vanishes when 9 — > tt/2. Even when 9 is very close to n/2, there- 
fore, the instability could come only from the effect of the across-slot gravity 
component. This effect is known to result in such a long- wavelength instability 



The viscous manifestation of linear steady instability in a vertical-slot geome- 
try seems to be limited to a more favorable background ratio between the op- 
posing component gradients. Indeed, enhancement of the relative background 
solute scale with respect to that in a LHSSS (Ras > 0, ji = 0) [21], taking 
place as /i increases from 0, is expected to be unfavorable for development of 
the LHSSS linear steady instability. In particular, the relative increase of the 
solute perturbation amplitude, associated with such solute scale enhancement, 
would result in a decrease of the horizontal density differences between respec- 
tive streamline points of the perturbation pattern illustrated in Fig. 2(a) of 
[22] . This was confirmed in trial computations for Ras = 30000 and A = 2. 

Essentially the same findings are reported for Ras = as well (Fig. ITTT) . when 
/j, increases from ~ 0.92 to ~ 0.98. The streamline perturbation pattern in Fig. 
HTIfa). which is similar to that in Fig. 2(a) of [22], could not be maintained as 
the relative solute perturbation amplitude grows when /i is increased towards 
1. Such pattern is thus found to undergo changes [Figs. ITlTb) and [TlTc)]. 
The changed streamline pattern in Fig. ITiTc), however, only enhances the 
unfavorable role of dissipation, which results in vanishing of the instability. 

The reported analysis of linear steady instability for Ras = was possible 
in view of the existence of a narrow interval of /i (~ 0.86 < fi <~ 0.98) for 
which A = 2 is unstable to steady infinitesimal disturbances. This interval has 
been uncovered in the present study. When /x is decreased below ~ 0.86, the 
linear steady instability disappears. Vertical shear motion, intensifying with 
decreasing /j, becomes an important factor in the nature of the background 
flow. This prevents the small-amplitude manifestation of differential gradi- 
ent diffusion. As fi exceeds ~ 0.98, on the other hand, the ratio between the 
background component scales becomes insufficiently favorable for the mani- 
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festation of small-amplitude steady convection. This ratio matters so much 
because of the effect of dissipation on the instability feedback. 

The quasisteady nature of oscillatory instability in the vertical-slot geometry 
suggests that the relative values of the component perturbations are critically 
important for the manifestation of such instability as well. In particular, the 
streamline patterns in Figs. [3(c) 13(g), and [3(h) are reminiscent of those in 
Fig. dU The maximal absolute values of solute perturbation in Figs. [3(c) 13(g), 
and [3(h) are however smaller, compared to the respective temperature per- 
turbation values, than those in Fig. [Hj Counter-propagation of the streamline 
disturbance constituents prevents such disturbance from coherently maintain- 
ing the maximal absolute values of the solute perturbation exceeding those of 
the temperature perturbation. The dominant part of either traveling stream- 
line perturbation pattern also interacts only with one of the sidewall boundary 
conditions. In contrast to this, the steady disturbance interacts with both such 
conditions. 



4-1.2 Effect of the across-slot gravity. 

When Ra > and 9 = it, different boundary conditions result in steady 
linear instability. The mechanism of this instability [2"0T2"T] has to be entirely 
responsible for small-amplitude steady convection for 9 G (ir/2, it) as well (Sec. 

At any 9 G (7r/2,7r) illustrated in Figs. 191 and [TUl the instability is restricted 
to the wave numbers that are smaller than a certain critical value. From the 
qualitative point of view, this is consistent with enhancement of the overall 
dissipation by the wave number increase. However, such largest unstable wave 
number decreases to zero as 9 approaches n/2. This does not seem to be a 
necessary consequence only of the vanishing of the across-slot gravity compo- 
nent. In addition, the increases of Ra c (k) in the vicinity of the largest unstable 
k for a given 9 in Figs. [9] and [10] are noticeably more abrupt than one could 
expect from the behavior of Ra c (k) at the same k and larger 9. The behavior 
of Ra c {k) is thus affected by an additional stabilizing factor. 

The curves in Figs. [9] and [10] are affected by the along-slot gravity component. 
In the presence of this gravity component, the steady perturbation maintained 
by the across-slot component would have to change its pattern. However, the 
along-slot gravity component alone fails to give rise to linear steady instability. 
Such perturbation change could thus be only stabilizing, for the resulting 
pattern would no longer be optimal with respect to the effect of the across- 
slot gravity component. In particular, the perturbation change could result 
in the across-slot cell path acquiring a slope with respect to the across-slot 
coordinate. This would enhance the dissipation of across-slot motion. Such 
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stabilizing effect is also most pronounced for small wavelengths, and its relative 
significance decreases as the wavelength increases. 

Another feature of Figs. M and [TU] is that the interval of unstable wave num- 
bers is longer for no-slip boundary conditions (Fig. [TU]) when lAn/2 < 9 < it. 
For 7r/2 < 9 < 1.3n/2, it is the stress-free conditions (Fig. [9]) that give rise 
to a longer interval of unstable wave numbers. In the latter case, all unstable 
wavelengths are large enough for the dissipation of along-slot perturbation 
motion to dominate the overall dissipation. The stress-free boundary condi- 
tions thus allow to maintain a longer interval of unstable wavelengths. As the 
slot orientation is far enough from vertical, however, the shortest unstable 
wavelengths are such that the dissipation of across-slot motion becomes rele- 
vant. The respective critical component gradients are then larger for no-slip 
boundary conditions. This factor proves to be more important for maintaining 
the instability of the shortest wavelengths than the dissipation decrease due 
to the stress-free boundary conditions. 

To obtain the dependence of Ra c (Q) on 9, 

Ra c {k) = Ra + Ra 2 k 2 + 

u(k) = k 2 (u + u 2 k 2 + ...), 

i(k) = t + t 2 k 2 + 

s(k) = s + s 2 k 2 + ... (32) 

were introduced into Eqs. ffT9l) — fl2T]) . t = and §o = const were thus obtained 
from the k° order of Eqs. f[2"0l and f[2~Tl) and the boundary conditions for to 
and s . Such t and s were then introduced into the k l order of Eq. f[T§j) . 
Along with either (stress-free) boundary conditions ffTT]) or (no-slip) boundary 
conditions (fT2"]) . this resulted in the respective expressions for uq. With the 
use of these expressions, the k 2 order of Eq. (|2TT) and the respective boundary 
conditions for s 2 lead to 

Ra c (0) = Ra = -Ra w J cos 9 = Ra*/ sin(x7r/2). (33) 

Here Ral is the value Ra c (0) would take on in the horizontal slot (9 = it): 
Ral = 120 for stress-free boundary conditions and Ra^ = 720 for the no- 
slip conditions [32]. Such values were first encountered in Rayleigh — Benard 
convection with the flux unstable gradient |33j . 

The derivation of Eq. (!33l) shows that the terms associated with the along- 
slot gravity component in the right-hand side of Eq. (fl9j) have no impact on 
Eq. (1331) . This is consistent with the above interpretations that the effect of 
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this gravity component is not among the causes of linear steady instability for 
9 G (7r/2, 7r) and that any role of this effect vanishes as k — > 0. Eq. (133|) was 
also found to accurately coincide with the numerical data underlying Figs. 
andOTOJ 

4-2 Finite- amplitude steady convection 
4-2.1 General. 

The specification of solute concentration at the middle points of the across- 
slot boundaries was found to numerically permit two pairs of periodic steady 
convective solutions arising from the conduction state. Within the selected 
along-slot period, one of such solution pairs is characterized at the onset of 
linear instability by two full counter-rotating cells. Except for the horizontal- 
slot geometry, the corresponding finite-amplitude solutions are generally re- 
flectionally asymmetric. This is the result of an asymmetry between their 
counter-rotating cells (see Sec. 14.2.21 below). 

In the other solution pair, a clockwise- (counterclockwise-)rotating cell in the 
center of the period interval is combined with parts of counter clockwise- 
(clockwise-)rotating cells near the across-slot boundaries. This pair possesses 
the reflection symmetry identified by LS[ip(— x, —y),t(—x, —y), s(—x, — y)} T = 
[ip(x, y), — t(x, y), — s(x, y)] T for all 9. The eigenvectors in Fig. [TTJ, in partic- 
ular, correspond to the solutions that are LS-symmetric for the illustrated 
along-slot interval. 

The only essential difference between the branch pairs is however the along-slot 
phase shift between their solutions. In particular, when linear steady instability 
is present for a given 9 (say for 9 = it or for 9 = as fi < 1), these solution pairs 
bifurcate from the conduction state at the Rayleigh numbers the difference 
between which is practically negligible. When the bifurcation is subcritical, 
the limit points are also formed at practically the same Ra. These pairs of 
branches can thus be viewed as standing for a single physically relevant pair 
of solution branches arising from the respective conduction state. 

The two pairs of convective solution branches were however found to un- 
dergo dissimilar secondary steady bifurcation phenomena, as Ra or 9 were 
varied. Such phenomena were detected by changes in the solution Jacobian 
sign, monitored during the continuation procedure. This dissimilarity could 
arise, in particular, from the different symmetries of the respective convective 
solutions, only one pair of which is LS'-symmetric. In terms of the flow struc- 
ture, however, the encountered secondary bifurcation phenomena were often 
found to give rise only to along-slot phase shifts of the respective convective 
patterns. The presence or absence of secondary steady bifurcations, and of the 
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instabilities associated with them, in real fluid could not thus be established 
with certainty in the framework of the present numerical approach. 

In addition, the stability of finite-amplitude steady convective branches to 
oscillatory disturbances was found to depend on which of the two numeri- 
cally present pairs of solutions is examined. This proved to hold even for the 
horizontal- slot geometry, where the flow pair with full counter-rotating cells 
also has a reflection symmetry. 

The prescribed across-slot-boundary values of s can prevent traveling-wave 
disturbances from being detected in the present formulation. However, finite- 
amplitude traveling-wave branches could intersect steady convective branches 
and thus affect the stability of the latter. Such findings have been previously 
reported for binary fluid convection [33] and for conventional double-diffusive 
convection [35]. The presence of traveling- wave branches also affects the sta- 
bility of reflectionally symmetric oscillatory branches [31] • This could have 
an additional effect on the stability of finite-amplitude steady solutions. The 
stability of steady convective branches to secondary oscillatory disturbances 
may therefore anyway change if the existence of traveling waves is properly 
allowed for. 

The solution structures reported below thus refer to the two pairs of steady 
convective branches arising in such numerical computations as one, as has been 
done for the periodic solution computations in [22|24II27] as well. The focus of 
this study is also on the physical mechanisms of conduction-state instability. 
For the reasons just mentioned, therefore, the issues of secondary (steady and 
oscillatory) instability of the computed finite-amplitude steady flows were left 
outside the scope of this work. 

4.2.2 0e [71-/2,71-]. 

For a large enough perturbation amplitude at 6 = tt/2, an appropriate stream- 
line pattern is maintained by the combined action of differential gradient dif- 
fusion and along-slot gravity [23]. This is possible despite the effects of dissi- 
pation, due to the energy coming from the nonlinear terms in Eqs. (pQ) — (J3J). 

As 9 changes from 7r to n/2, the mechanism of differential gradient diffu- 
sion discussed in [201121] thus gradually gives way to that considered in |24J. 
Via supercritical bifurcation, the former mechanism gives rise to a convective 
steady pattern of symmetric counter-rotating cells. The steady convection pat- 
tern resulting from the finite-amplitude mechanism in [23] has however only 
clockwise-rotating cells. The along-slot gravity component has an asymmetric 
effect on the senses of rotation of the counter-rotating cells. 

In particular, let a perturbation pattern of symmetric counter-rotating cells 
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be generated in a vertical slot with two background horizontal gradients and 
the higher component values at the left sidewalk The clockwise-rotating cells 
would then create locally unstable solute perturbation stratifications. Such 
stratifications would however be stable where the cells rotate counterclock- 
wise. Due to differential diffusion (whether of gradient or classical nature), 
the clockwise-rotating cells have to intensify with respect to the cells rotating 
counterclockwise [22,23,24,26]. The clockwise-rotating cells are thus dominant 
in the finite-amplitude steady pattern. 

At the onset of linear steady instability, the counter-rotating cells are bound 
to be symmetric. The asymmetry between them could thus arise, as 9 changes 
from 7i to 7r/2, only with a finite convection amplitude. The latter is expected 
to arise when the Rayleigh number increases above the linear instability onset. 
The effect of the along-slot gravity component is itself also enhanced as the 
Rayleigh number is increased. For Ra = 31000 and A = 2, in particular, 
the flow patterns are thus entirely dominated by the clockwise-rotating cells 
when 9 is still greater than 1.97r/2 [Figs. [T27 a) — H2T d)]. The remnants of the 
counterclockwise-rotating cells also vanish for 9 G (1.77r/2, 1.757r/2) [Figs. 

For 9 G [7r/2,7r], differential gradient diffusion underlies the convective mech- 
anisms resulting from both the along-slot and the across-slot gravity compo- 
nents. The mechanism at 9 = tt/2, however, leads only to nonlinear steady 
instability. As 9 changes from tt to tt/2, primary control over finite-amplitude 
steady convection is thus continuously transformed from the convective mech- 
anism arising from the across-slot gravity component to that associated with 
the along-slot component. Such linear instability can however be due to the 
decreasing across-slot gravity component alone. At some 9[ « (1.6 — 1.7)7r/2 
for A = 2], this renders the steady bifurcation subcritical, as seen from Figs. 
rrSTa) andflUVb). As the across-slot gravity component further decreases with 
9 being decreased [to 9 G (1.37r/2, lAn/2) for A = 2], the linear steady insta- 
bility vanishes [Figs. [TSTc) and [TBT d)]. 



4.2.3 0e [0,7r/2]. 

For 9 G [0,7r/2], linear steady instability can arise for // = 1 neither from 
the effect of the along-slot gravity component nor from that of the across- 
slot component. The latter effect favors amplitude growth only of the small- 
amplitude convective cells whose sense of rotation changes periodically in time 
(Sees. [3T2TTl and [3~.3.1|) . In addition, differential gradient diffusion is not among 
the causes of finite-amplitude steady convection at 9 = 0. Rather, it plays only 
a stabilizing role in the mechanism of such convection [27]. 

The nature of finite-amplitude steady convective flows at 9 = is thus ba- 
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sically different from that at 9 = n/2. As the convection amplitude becomes 
large enough, with the Rayleigh number being increased, such dissimilarity 
in the nature of convection prevents continuous transformation of one of the 
convective steady flows into the other. In particular, a region of hysteresis in 
9 is formed for Ra >~ 20000 (Fig. [Til). 

As 9 increases from 0, finite-amplitude steady convection at Ra = 31000 re- 
mains under primary control of the effect of the across-slot gravity component. 
Playing a stabilizing role in this effect [27], differential gradient diffusion re- 
duces the relative significance of its destabilizing role in the emerging effect of 
the along-slot gravity component. The counterclockwise-rotating cells thus re- 
tain their horizontal-slot pattern along the whole branch A2h [Figs. [ToT a) and 
H51(b)]. For Ra = 31000 and A = 2, this branch extends to about 9 « 0.51tt/2 
[Fig. He)]. 

With 9 decreasing from tt/2, finite- amplitude steady convection remains un- 
der primary control of the effect of the along-slot gravity component. This 
effect is underlain by differential gradient diffusion [2"2~|T2"1] . Such process re- 
duces the effectiveness of the convective mechanism by means of which the 
emerging across-slot gravity component acts [27]. This gravity component 
thus fails to have a significant effect on the convective pattern (i.e., to gen- 
erate counterclockwise- rotating cells). The whole branch A2 V , extending to 
9 pa 0.387r/2 for Ra = 31000 and A = 2, is therefore characterized only by 
clockwise-rotating cells [Fig. [ToT d)]. Along branch A2 hv [Fig. [TSlfc)]. such con- 
vection pattern is transformed into that with counter-rotating cells [Figs.fToTb) 
and[15Kc)]. 



5 Summary and concluding remarks 

This work provides a comprehensive understanding of the manifestation of 
double-component convection due to the effects of different component bound- 
ary conditions in an infinite slot having diverse orientations to the gravity. To 
eliminate other physical effects, the primary focus has been on the problem 
with compensating gradients of the components and equal component diffu- 
sivities. Small-amplitude oscillatory convection as well as small- and finite- 
amplitude steady convection have been considered. Although the study was 
primarily focused on 2D convection, small-amplitude 3D effects were also ex- 
amined where the obtained physical insight suggested their potential relevance. 

The main aspect of the provided understanding is a consistent interpretation 
of the physics underlying small- amplitude oscillatory convection. For all slot 
orientations, the physical nature of such convection both in inviscid and in 
viscous fluid proved to be underlain by differential gradient diffusion. However, 
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the specific mechanisms by means of which the manifestation of this effect 
takes place vary significantly with the slot orientation. Their features also 
depend substantially on whether the involved fluid is of inviscid or viscous 
nature. 

In the horizontal slot at 9 = with inviscid fluid and an independently pre- 
scribed flux stratification, Ra s , the marginal-stability boundaries behave con- 
sistently with the physical mechanism of oscillatory instability first highlighted 
in [T5]. In particular, the most unstable wave number is zero and the instabil- 
ity arises within a finite wave number interval. The oscillation frequency also 
decreases with decreasing wave number and becomes zero at k = 0. By means 
of the long-wavelength expansion, the critical fixed-value stratification, Ra c , 
and the group velocity, du c /dk, for k = were expressed in terms of Ra s . 

For the compensating background gradients in inviscid fluid, the most unsta- 
ble wave number remains equal to zero up until the slot orientation becomes 
nearly vertical (9 = 7r/2). The analytical relation between Ra c and dui c /dk at 
k = was derived for such inclined slot at any 9 £ [0, 7r/2) as well. For 9 grow- 
ing from to 7r/2, the oscillatory instability is however increasingly affected by 
the along-slot gravity component. The effect of this gravity component is asso- 
ciated with the mechanism of steady instability identified in [22.24J for viscous 
fluid. The oscillatory marginal-stability curves for 9 £ (0, 7r/2) were found to 
behave consistently with the assumption that the corresponding mechanism in 
inviscid fluid is effective at any Ra ^ 0. In particular, the oscillatory marginal- 
stability boundary was found to transform into the steady one with Ra c (k) = 
as 9 -> tt/2. 

In the horizontal slot (9 = 0) with viscous fluid, the vicinity of zero wave 
number is no longer the region of most unstable k. u> c (k) also tends to a nonzero 
value as k — > 0. These findings are still consistent with the interpretation of 
the oscillatory instability for inviscid fluid. In particular, the dissipation of 
along-slot perturbation motion was found to be part of the feedback in the 
horizontal-slot oscillatory instability mechanism. Increase of the instability 
wavelength thus reduces the efficiency of such feedback. As k approaches 0, 
this stabilizing factor becomes more important than the destabilizing decrease 
of the overall dissipation. 

The suggested interpretation of the role played by the dissipation of along-slot 
motion in the horizontal-slot oscillatory instability feedback apparently applies 
to the classical double-diffusion as well. In particular, zero wave number has 
to be most unstable in the diffusive regime of conventional double-diffusive 
convection when both background gradients in inviscid fluid are maintained 
by flux boundary conditions. Indeed, Fig. 4(b) in [36] seems to suggest that 
the most unstable wavelength tends to infinity and the oscillation frequency 
tends to zero as Pr — > 0. For a finite Pr, however, both the most unstable 
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wavelength and the corresponding oscillation frequency in Fig. 4(b) of [36] are 
finite. 

For relatively large X y , most unstable disturbances of the oscillatory instability 
in an inclined slot [0 6 (0, 7r/2)], where the across-slot gravity component is 
involved in driving the instability, are of a 3D nature. This was also suggested 
to result from the effect of dissipation of along-slot perturbation motion. In 
particular, the feedback efficiency increases when A decreases due to orien- 
tation change of the 2D axis of rotation of a perturbation cell. All relevant 
findings proved to be consistent with such interpretation of the 3D instability 
nature. For small 9 > 0, this interpretation was also found to be a manifes- 
tation of more general conditions for three-dimensionality of instability. Such 
conditions were formulated independently of the physics underlying a specific 
problem. 

Another important outcome of this study is identification of the oscillatory in- 
stability mechanism arising from the effect of the along-slot gravity component 
in viscous fluid. The reflectionally symmetric perturbation mode of this insta- 
bility consists of two counter-propagating patterns either of which is dominant 
near one of the slot sidewalls. The counter-traveling convective streamline pat- 
terns give rise to such feedback distribution of the component perturbations 
as generates their amplitude growth and a nearly steady density difference 
across the slot. This density difference sustains the counter-propagation of the 
convective patterns. Such growth of the convection amplitude and formation 
of the across-slot density difference are underlain by differential gradient dif- 
fusion. In a quasisteady sense, the manifestation of such differential diffusion 
is analogous to that in the respective steady instability [22~f24"] . 

The identified mechanism of vertical-slot oscillatory instability is apparently 
also relevant under the conditions other than those considered in this work. In 
particular, Fig. 5 in [37J suggests that oscillatory instability arises in the verti- 
cal slot filled with a porous medium even when different boundary conditions 
could be the only cause of this instability. In addition, the instability mech- 
anism identified in the present study may also be effective when the role of 
differential gradient diffusion in it is played by the classical double-diffusion. 
Indeed, the directions in which counter-propagating flow patterns travel in 
Fig. 8 of [37] seem to be consistent with the respective directions in Fig. [7J of 
the present work. The eigenfunctions of traveling waves illustrated in Fig. 7 of 
[37] are also dominated by their patterns near the respective slot boundaries. 

The quasisteady nature of oscillatory instability in the vertical slot suggests 
that the dissipation of along-slot perturbation motion is part of the instability 
feedback, as also highlighted for the respective mechanism of steady instability. 
Since such a dissipation is part of the oscillatory instability feedback at 6 = 
as well, the same shape of the oscillatory marginal-stability boundaries for all 
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9 G [0, 7r/2] is thus explained. In contrast to small- amplitude steady vertical- 
slot convection, the manifestation of such quasisteady oscillatory convection 
is permitted due to a relatively more favorable ratio between the maximal 
absolute values of the component perturbations. 

As 9 exceeds tt/2, the mechanism of linear steady instability arising from 
the effect of the across-slot gravity component [20||2T] becomes increasingly 
involved in the manifestation of oscillatory instability due to the along-slot 
component. Behavior of the oscillatory marginal-stability curves above 9 « 
1.57r/2 was thus found to suggest their transformation into the respective 
curves characterized by a zero real eigenvalue. 

The effect of the across-slot gravity component identified in [201121] is the only 
cause of linear steady instability for 9 G (7r/2, tt]. Driven by differential gradi- 
ent diffusion, such instability arises within a finite interval of wave numbers 
that extends from the most unstable zero wave number. Obtained from the 
long-wavelength expansion, Ra c at k = is not affected by the along-slot 
gravity component at all. However, the along-slot gravity component still has 
a stabilizing effect on small-amplitude steady perturbations with k > 0. This 
effect is manifested in that the interval of unstable wave numbers shrinks to 
the vicinity of k = and vanishes with 9 decreasing from n to tt/2. 

The effect of the along-slot gravity component still gives rise to finite-amplitude 
steady convection. In contrast to the finite-amplitude pattern of symmetric 
counter-rotating convective cells at 9 = n, the convection pattern at 9 = it/ 2 
comprises only clockwise-rotating cells. Underlain by differential gradient dif- 
fusion, the basic nature of finite-amplitude convective steady flows at 9 = tt/2 
is however the same as that at 9 = tt. Continuous transformation between 
these flows was thus demonstrated to take place at least up to Ra = 31000. 

In contrast to 9 G [tt/2,tt], finite-amplitude convective steady flows at 9 = 
could not be continuously transformed into those at 9 = n/2 for Ra = 31000. 
Rather, a region of hysteresis in 9 was demonstrated to develop above Ra ~ 
20000. Such hysteresis phenomenon is a manifestation of the basic nature 
of finite-amplitude steady convection at 9 = being different from that at 
9 = tt/2. In particular, differential gradient diffusion is a stabilizing factor with 
respect to the effect of the across-slot gravity component on finite-amplitude 
steady convection for 9 G [0, tt/2). It is the opposite roles of this physical 
process in the effects of the across-slot and along-slot gravity components 
that prevents continuous transformation between the finite-amplitude steady 
flows at 9 = and 9 = tt/2. 

The reported mechanisms of double- component convection due to different 
component boundary conditions thus comprise two basic ingredients. These 
are the disparity between component stratifications caused by a convective 
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perturbation and the resulting feedback that maintains amplitude growth of 
such a perturbation. For all slot orientations, the feedback in small-amplitude 
oscillatory and steady convection is underlain by differential gradient diffusion. 
In particular, such universality of the nature of small-amplitude convection is 
manifested in the transformations of oscillatory marginal-stability curves into 
steady ones. Finite-amplitude steady convection can however arise due to two 
dissimilar feedbacks. The role of differential gradient diffusion is destabilizing 
in one of them and stabilizing in the other. 

The disparity between component gradients in perturbed state is thus the only 
universal feature of the effects of different boundary conditions. This feature 
is however also independent of the physical nature of the forces responsible 
for the feedback. It could therefore be expected to give rise to instabilities 
when the forces generated by components are other than the buoyancy forces. 
Analysis of such possibilities could provide a further insight into the effects of 
different boundary conditions. 
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tf^ = 5S-5T 
t = -Ra/2 Sv=-Ra s 



t = Ra/2 Sv=-Ra s 



Fig. 1. Geometry of the problem. 5p = p — po is the variation of the (dimensionless) 
density, p, due to the variations SS and 5T of solute concentration S and tempera- 
ture T = (Ti + T2)/2 + t with respect to their reference values, at which the density 
is po; T\ and T2 are the boundary temperatures, a = gd 3 /KV. Pr = v/k = 6.7, 
Le = kt I us = 1; k>t and ( = k) are the component diffusivities. The fluid is of 
the Boussinesq type. 
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Fig. 2. 6 = 0. Inviscid fluid. Curves of the marginal linear stability to 2D oscil- 
latory disturbances for independently prescribed Ra s , Ra c (k) and u) c (k); Le = 1. 
(a) the solid lines: Ra s = 50000, the dashed lines: Ra s = 20000, the dotted lines: 
Ra s = 10000, the dash-dot lines: Ra s = 5000; (b) the solid lines: Ra s = 1000, the 
dashed lines: Ra s = 500, the dotted lines: Ra s = 200, the dash-dot lines: Ra s = 50. 
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(a) 



(b) 



(c) 



Fig. 3. Inviscid fluid. Curves of the marginal linear stability to 2D (k z = 0) oscil- 
latory disturbances, Ra c (k) and u> c (k), and to such 3D (k z > 0) disturbances with 
different values of k z , Ra c {k y ) and uj c {k y ); fi = 1, Le = 1. The numbers in (b) — (e), 
where the disturbances with k z > were also considered, give the values of k z cor- 
responding to their nearest curve they do not intersect, (a) 6 = 0; (b) 6 = ir/8; 
(c) 6 = vr/4; (d) 6 = 3vr/8; (e) 6 = 0.95vr/2; (f) the dotted lines: = 0.99vr/2, the 
dashed lines: 9 = 0.995vr/2, the solid lines: 9 = 0.999vr/2. 
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(d) 



(e) 



(f) 



Fig. 3. Inviscid fluid. Curves of the marginal linear stability to 2D {k z = 0) oscil- 
latory disturbances, Ra c (k) and ui c {k), and to such 3D {k z > 0) disturbances with 
different values of k z , Ra c (k y ) and to c (k y ); fx = 1, Le = 1. The numbers in (b) — (e), 
where the disturbances with k z > were also considered, give the values of k z cor- 
responding to their nearest curve they do not intersect, (a) 9 = 0; (b) 9 = ir/8; 
(c) 6 = vr/4; (d) 6 = 3vr/8; (e) 9 = 0.95vr/2; (f) the dotted lines: 9 = 0.99vr/2, the 
dashed lines: 9 = 0.995vr/2, the solid lines: 9 = 0.999vr/2. 
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(a) 



(b) 



(c) 



Fig. 4. Viscous fluid and stress-free slot boundaries. Curves of the marginal linear 
stability to 2D [k z = 0) oscillatory disturbances, Ra c (k) and u c (k), and to such 
3D (k z > 0) disturbances with different values of k z , Ra c (k y ) and LO c (k y ); ji = 1, 
Pr = 6.7, Le = 1. The numbers in (b) — (e), where the disturbances with k z > 
were also considered, give the values of k z corresponding to their nearest curve they 
do not intersect, (a) 6 = 0; (b) 6 = vr/8; (c) 6 = vr/4; (d) 6 = 3vr/8; (e) 9 = vr/2; 
(f) the dash-dot lines: 6 = 57r/8, the dotted lines: 9 = 1.357r/2, the dashed lines: 
9 = 1.45vr/2, the solid lines: 9 = 1.49vr/2. 
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Fig. 4. Viscous fluid and stress-free slot boundaries. Curves of the marginal linear 
stability to 2D {k z = 0) oscillatory disturbances, Ra c (k) and co c (k), and to such 
3D (k z > 0) disturbances with different values of k z , Ra c (k y ) and LO c (k y ); ji = 1, 
Pr = 6.7, Le = 1. The numbers in (b) — (e), where the disturbances with k z > 
were also considered, give the values of k z corresponding to their nearest curve they 
do not intersect, (a) 9 = 0; (b) 9 = vr/8; (c) 9 = it/A; (d) 6 = 3vr/8; (e) 6 = vr/2; 
(f) the dash-dot lines: 6 = 5ir/8, the dotted lines: 9 = 1.357r/2, the dashed lines: 
9 = 1.45tt/2, the solid lines: 9 = 1.49vr/2. 
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Fig. 5. Viscous fluid and no-slip slot boundaries. Curves of the marginal linear 
stability to 2D {k z = 0) oscillatory disturbances, Ra c (k) and co c (k), and to such 
3D (k z > 0) disturbances with different values of k z , Ra c (k y ) and LO c (k y ); [i = 1, 
Pr = 6.7, Le = 1. The numbers in (b) — (e), where the disturbances with k z > 
were also considered, give the values of k z corresponding to their nearest curve they 
do not intersect, (a) 9 = 0; (b) 9 = vr/8; (c) 9 = it/A; (d) 9 = 3vr/8; (e) 9 = vr/2; 
(f) the dash-dot lines: 9 = 5ir/8, the dotted lines: 9 = 1.3n/2, the dashed lines: 
9 = 1.47T/2, the solid lines: 9 = 1.48vr/2. 
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(e) 



(f) 



Fig. 5. Viscous fluid and no-slip slot boundaries. Curves of the marginal linear 
stability to 2D {k z = 0) oscillatory disturbances, Ra c (k) and co c (k), and to such 
3D (k z > 0) disturbances with different values of k z , Ra c (k y ) and LO c (k y ); ji = 1, 
Pr = 6.7, Le = 1. The numbers in (b) — (e), where the disturbances with k z > 
were also considered, give the values of k z corresponding to their nearest curve they 
do not intersect, (a) 9 = 0; (b) 9 = vr/8; (c) 9 = it/A; (d) 9 = 3vr/8; (e) 9 = vr/2; 
(f) the dash-dot lines: 9 = 5ir/8, the dotted lines: 9 = I.Stt/2, the dashed lines: 
9 = 1.47T/2, the solid lines: 9 = 1.48vr/2. 

46 




Fig. 6. 9 = 0. Viscous fluid and no-slip slot boundaries. Perturbation temporal be- 
havior throughout the oscillation period t p ~ 2A5t (5t = 0.05) just beyond the 
onset of 2D oscillatory instability. It was obtained from the numerical simulation 
of evolution of the linearized Eqs. ([I]) — (Ras = 0) in response to the initial 
disturbance proportional to the background state after initial time Tj ~ 19000 has 
passed; A = 2, [i = 1, Ra = 30060, Pr = 6.7, Le = 1. With this Tj, all pertur- 
bation modes other than the unstable mode (r p ~ 24<5r) are practically negligible. 
ip: perturbation streamlines; s: isolines of solute concentration perturbation; t: per- 
turbation isotherms. The actual relative values of the streamfunction perturbation 
are equal to 10 -3 times the respective values in the figure. The solid and dashed 
streamlines designate the clockwise and counterclockwise rotation and are equally 
spaced within the positive and negative streamfunction intervals, respectively. The 
solid and dash-dot isolines of the component perturbations designate the positive 
and negative component perturbation intervals, respectively, (a) r = Tj + 3<5t; (b) 
T = n + 65t; (c) T = n + 85t; (d) r = n + 10«5r; (e) r = n + 12<5r; (f ) r = n + Udr; 
(g) T = n + 17Jr; (h) r = Ti + 205r; (i) r = n + 22<5r; (j) r = n + 24Jr. 
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Fig. 6. 9 = 0. Viscous fluid and no-slip slot boundaries. Perturbation temporal be- 
havior throughout the oscillation period t p ~ 2A5t (5t = 0.05) just beyond the 
onset of 2D oscillatory instability. It was obtained from the numerical simulation 
of evolution of the linearized Eqs. ([I]) — (Ras = 0) in response to the initial 
disturbance proportional to the background state after initial time Tj ~ 19000 has 
passed; A = 2, [i = 1, Ra = 30060, Pr = 6.7, Le = 1. With this Tj, all pertur- 
bation modes other than the unstable mode (r p ~ 24<5r) are practically negligible. 
ip: perturbation streamlines; s: isolines of solute concentration perturbation; t: per- 
turbation isotherms. The actual relative values of the streamfunction perturbation 
are equal to 10 -3 times the respective values in the figure. The solid and dashed 
streamlines designate the clockwise and counterclockwise rotation and are equally 
spaced within the positive and negative streamfunction intervals, respectively. The 
solid and dash-dot isolines of the component perturbations designate the positive 
and negative component perturbation intervals, respectively, (a) r = Tj + 3<5t; (b) 
T = n + 65t; (c) T = n + 85t; (d) r = n + 10«5r; (e) r = n + 12<5r; (f ) r = n + Udr; 
(g) T = n + 17Jr; (h) r = Ti + 205r; (i) r = n + 22<5r; (j) r = n + 24Jr. 
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Fig. 7.8 = 7r/2. Viscous fluid and no-slip slot boundaries. Perturbation temporal 
behavior throughout the oscillation period t p ~ 345t (<5t = 0.05) just beyond the 
onset of 2D oscillatory instability. It was obtained from the numerical simulation of 
evolution of the linearized Eqs. ([I]) — (j4|) (Ras = 0) in response to the initial distur- 
bance proportional to the background state after initial time Tj ~ 2000 has passed; 
A = 2, n = 1, Ra = 17411, Pr = 6.7, Le = 1. With this Tj, all perturbation modes 
other than the unstable mode (t p ~ 345t) are practically negligible, tp: perturba- 
tion streamlines; s: isolines of solute concentration perturbation; t: perturbation 
isotherms. The actual relative values of the streamfunction perturbation are equal 
to 10 -3 times the respective values in the figure. The solid and dashed streamlines 
designate the clockwise and counterclockwise rotation and are equally spaced within 
the positive and negative streamfunction intervals, respectively. The solid and dash- 
dot isolines of the component perturbations are equally spaced within the positive 
and negative component perturbation intervals, respectively, (a) r = Tj + 35t; (b) 
r = Ti + 75r; (c) r = n + lUr; (d) r = n + lbdr; (e) r = Tj + 19<5r; (f) r = Tj + 225r; 
(g) T = Ti + 25<5r; (h) r = n + 285t; (i) T = n + 31<5r; (j) r = n + 345t. 
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Fig. 7.8 = ir/2. Viscous fluid and no-slip slot boundaries. Perturbation temporal 
behavior throughout the oscillation period t p ~ 345t (<5t = 0.05) just beyond the 
onset of 2D oscillatory instability. It was obtained from the numerical simulation of 
evolution of the linearized Eqs. ([I]) — (j4|) (Rag = 0) in response to the initial distur- 
bance proportional to the background state after initial time Tj ~ 2000 has passed; 
A = 2, n = 1, Ra = 17411, Pr = 6.7, Le = 1. With this Tj, all perturbation modes 
other than the unstable mode (t p ~ 345t) are practically negligible, tp: perturba- 
tion streamlines; s: isolines of solute concentration perturbation; t: perturbation 
isotherms. The actual relative values of the streamfunction perturbation are equal 
to 10 -3 times the respective values in the figure. The solid and dashed streamlines 
designate the clockwise and counterclockwise rotation and are equally spaced within 
the positive and negative streamfunction intervals, respectively. The solid and dash- 
dot isolines of the component perturbations are equally spaced within the positive 
and negative component perturbation intervals, respectively, (a) r = Tj + 35t; (b) 
r = Ti + 75r; (c) r = n + lUr; (d) r = n + lbdr; (e) r = Tj + 19<5r; (f) r = Tj + 225r; 
(g) T = Ti + 25<5r; (h) r = n + 285t; (i) T = n + 31<5r; (j) r = n + 345t. 
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Fig. 8. The same as in Fig. [7J p: isolines of the perturbation in s — t. The solid 
and dash-dot isolines are equally spaced within the positive and negative density 
perturbation intervals, respectively. 
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Fig. 9. Viscous fluid and stress-free slot boundaries. Curves of the marginal linear 
stability to 2D steady disturbances, Ra c (k), for different 6 = (l + x)vr/2, < % < 1; 
fi = 1, Le = 1. The data for the wave numbers very close to k = (where stable 
data were numerically difficult to obtain) are not presented. 
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Fig. 10. Viscous fluid and no-slip slot boundaries. Curves of the marginal linear 
stability to 2D steady disturbances, Ra c (k), for different 8 = (l + x)vr/2, < % < 1; 
H = 1, Le = 1. The data for the wave numbers very close to k = (where stable 
data were numerically difficult to obtain) are not presented. 
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Fig. 11. 9 = tt/2. Viscous fluid and no-slip slot boundaries. Singular eigenvectors 
corresponding to the wavelength A = 2 at the onset of small-amplitude steady con- 
vection within a narrow interval of \x where the linear steady instability for A = 2 
arises; Pr = 6.7, Le = 1. ip: perturbation streamlines; s: isolines of solute concentra- 
tion perturbation; t: perturbation isotherms. The variables are nondimensionalized 
as in Eqs. ([T]) — (jU) (Ras = 0). The actual relative values of the streamfunction per- 
turbation are equal to 10 -3 times the respective values in the figure. The solid and 
dashed streamlines designate the clockwise and counterclockwise rotation and are 
equally spaced within the positive and negative streamfunction perturbation inter- 
vals, respectively. The solid and dash-dot isolines of the component perturbations 
designate the positive and negative component perturbation intervals, respectively, 
(a) fj, = 0.92, Ra c = 6139; (b) fj, = 0.96, Ra c = 8900; (c) fi = 0.98, Ra c = 22469. 
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Fig. 12. Viscous fluid and no-slip slot boundaries. Finite-amplitude convective 
steady flows representing the higher-amplitude branch, A2, for different 9 € (vr/2, 7r); 
A = 2, /x = 1, i?a = 31000, Pr = 6.7, Le = 1. The steady solutions were obtained 
from Eqs. (P) — (i?as = 0). The across- and along-slot coordinate axes are di- 
rected rightwards and upwards, respectively. ^: streamlines; S: isolines of solute 
concentration; t: isotherms. The actual values of t and s are equal to 10 4 times the 
respective values in the figure. The solid and dashed streamlines designate the clock- 
wise and counterclockwise rotation and are equally spaced within the positive and 
negative streamfunction intervals, respectively. The solid and dash-dot component 
isolines designate the positive and negative component values, respectively. The zero 
isotherms are designated by the dotted lines, (a) 9 = 1.997r/2; (b) 9 = 1.977r/2; (c) 
9 = 1.95tt/2; (d) 9 = 1.94tt/2; (e) 9 = 1.9vr/2; (f) 9 = 1.75vr/2; (g) 9 = 1.7tt/2; (h) 
9 = l.lvr/2. 
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Fig. 12. Viscous fluid and no-slip slot boundaries. Finite- amplitude convective 
steady flows representing the higher-amplitude branch, A2, for different 9 S (vr/2, 7r); 
A = 2, fi = 1, Ra = 31000, Pr = 6.7, Le = 1. The steady solutions were obtained 
from Eqs. ([I]) — (J4j) (Rag = 0). The across- and along-slot coordinate axes are di- 
rected rightwards and upwards, respectively. ^: streamlines; S: isolines of solute 
concentration; t: isotherms. The actual values of t and s are equal to 10 4 times the 
respective values in the figure. The solid and dashed streamlines designate the clock- 
wise and counterclockwise rotation and are equally spaced within the positive and 
negative streamfunction intervals, respectively. The solid and dash-dot component 
isolines designate the positive and negative component values, respectively. The zero 
isotherms are designated by the dotted lines, (a) 9 = 1.997r/2; (b) 9 = 1.977r/2; (c) 
9 = 1.95vr/2; (d) 9 = 1.94vr/2; (e) 9 = 1.9vr/2; (f) 6 = 1.75vr/2; (g) 6 = 1.7vr/2; (h) 
9 = 1.17T/2. 
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Fig. 13. Viscous fluid and no-slip slot boundaries. Schematic diagrams illustrat- 
ing structures of the steady flows representing onset of the linear and finite-am- 
plitude steady instability of the conduction state as Ra is increased for different 
9 = (x + l)vr/2, -1 < x < 1; A = 2, \i = 1, Pr = 6.7, Le = 1. $ is an abstract 
measure of the steady flows that distinguishes between different steady solutions, 
specifies the location of the singularities (limit points and symmetry-breaking bi- 
furcations), and represents the flows arising from a symmetry-breaking bifurcation 
as a single branch. The solid (dashed) lines represent the branches expected to be 
linearly stable (unstable) to the disturbances associated with the eigenvalues that 
give rise to the steady instability of the conduction state. Secondary oscillatory and 
steady bifurcations, if any, are not shown. B is the symmetry-breaking bifurcation 
standing for the steady linear stability boundary for wave number k = ir (A = 2). 
L is the limit point. Al and A2 are the lower- and higher-amplitude branches as- 
sociated with the limit point, respectively, (a) 9 = 1.87r/2; (b) 9 = 1.67r/2; (c) 
9 = 1.4vr/2; (d) 9 = 1.3rr/2; (e) 9 = 0.7vr/2; (f) 9 = 0.3vr/2. 
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Fig. 14. Viscous fluid and no-slip slot boundaries. Schematic diagrams illustrating 
structures of the steady flows of the higher-amplitude convective branch, A2, with 
the minimal along-slot period A = 2 as 9 = r\n j2 is varied for different Ra, rj G [0, 1]; 
fj, = 1, Pr = 6.7, Le = 1. $ is an abstract measure of the steady flows that 
distinguishes between different steady solutions and specifies the location of the 
singularities (limit points). The dashed lines imply that the respective branches are 
linearly unstable to steady disturbances. A2h and A2 V are the convective branches 
continuously transformed into the higher-amplitude branches in the horizontal and 
vertical slot, respectively. A2h v is the linearly unstable (to steady disturbances) 
branch connecting branches A2^ and A2 V via limit points LI and L2. Secondary 
instability, if any, of the depicted steady branches to the oscillatory and steady 
disturbances other than the (steady) disturbances destabilizing branch A2^ v is not 
shown, (a) Ra = 19000; (b) Ra = 22000; (c) Ra = 31000. 
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Fig. 15. Viscous fluid and no-slip slot boundaries. Finite- amplitude convective 
steady flows representing the higher-amplitude branch, A2, for different 9 £ (0, vr/2); 
A = 2, fj, = 1, .Ra = 31000, Pr = 6.7, Le = 1. The steady solutions were obtained 
from Eqs. (P) — (|3|) (Rag = 0). The across- and along-slot coordinate axes are di- 
rected rightwards and upwards, respectively. streamlines; S: isolines of solute 
concentration; t: isotherms. The actual values of t and s are equal to 10 4 times 
the respective values in the figure. The solid and dashed streamlines designate the 
clockwise and counterclockwise rotation and are equally spaced within the positive 
and negative streamfunction intervals, respectively. The solid and dash-dot com- 
ponent isolines designate the positive and negative component values, respectively. 
The zero isotherms are designated by the dotted lines, (a) 6 = 0.2-7r/2, branch A2^\ 
(b) 6 = 0.5vr/2, branch A2 hv ; (c) 6 = 0.4vr/2, branch A2 hv ; (d) 9 = 0.5vr/2, branch 
A2 V . 
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